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1.  Introduction 


The  overall  objectives  of  this  project  can  be  summarized  as  follows: 

(i)  To  obtain  experimental  justification  for  a  forward-scattering,  short  pulse  synthetic 
array  system  concept  towards  detection  and  possible  imaging  of  shallow  subsurface  targets  such 
as  anti-tank  and  anti-personnel  mines.  Consequently,  experimental  studies  (at  L-  and  S-  bands) 
constituted  the  primary  thrust  of  this  program  and  are  highlighted  in  the  body  of  the  report. 
Extensive  field  data  was  obtained  on  several  test  objects  as  well  as  mine  casings  supplied  by 
MERADCOM  when  buried  in  an  experimental  test-bed  (constructed  under  this  project)  under  a 
variety  of  conditions  at  the  Valley  Forge  Research  Center's  Experimental  facility,  operated  by 
the  Moore  School  of  Electrical  Enginnering,  University  of  Pennsylvania.  The  details  of  the 
experimental  work  are  outlined  in  [1-5].  The  promise  and  the  problems  of  the  experimental  set¬ 
up  employed  in  the  study  using  a  bistatic,  forward-scattering  configuration  is  described  in  the 
section  on  experimental  work. 

(ii)  Based  on  the  experimental  evidence  obtained,  develop  preliminary  propagation 
models  and  signal  processing  tools  for  location,  identification  and  classification  of  buried 
objects.  The  product  of  this  effort  include  sofware  routines  for  signal/image  processing  and 
image  display,  which  are  included  in  the  Appendix  along  with  examples  of  recent  images 
obtained  of  the  mine  casings.  Other  imaging  results  can  be  found  in  [1-5].  The  modest  success 
achieved  in  dry  soil  conditions  indicates  the  need  for  improved  modeling  of  wave  propagation 
and  interaction  at  the  interfaces,  to  be  used  by  subsequent  enhanced  signal  processing  tools  to 
obtain  high-resolution  imaging. 


2.  Experimental  Work 


2.1  Introduction.  General  System  Considerations 

The  goal  of  this  part  of  the  program  is  to  study  the  application  of  a  short  pulse,  high 
angular  resolution  array  system,  focused  in  the  near-field,  to  the  detection  and  imaging  of 
shallow  sub-surface  land  mines.  The  method  of  implementation  is  to  be  synthetic  aperture 
digital  beamforming,  in  which  the  array  element  outputs  are  sequentially  sampled,  and  the 
beamfarming,  signal  processing,  and  imaging  are  performed  in  the  system  computer. 

The  major  obstacles  to  achieving  these  goals  are  the  varying  dielectric  constant  and 
conductivity  of  typical  soil  types.  Land  mines  are  typically  only  a  few  centimeters  below  the 
soil  surface.  The  relatively  high  dielectric  constant  of  the  soil  will  cause  a  large  ground  surface 
reflection  which  may  be  larger  than  the  return  from  a  target,  and  both  returns  may  be  in  the  same 
range  bin;  even  a  pulse  as  short  as  one  nanosecond  has  a  range  resolution  in  sand  of  10  cm  for 
backscatter,  20  cm  for  forward  scatter,  and  probably  could  not  resolve  a  shallow  target  from  the 
surface  reflection.  The  finite  conductivity  of  the  soil  may  cause  a  large  attenuation  of  the  signal 
at  higher  frequencies,  as  well  as  phase  variations  in  the  signal  return.  Also,  target  dielectric 
constant  may  be  similar  to  that  of  the  soil  making  detection  difficult.  Finally,  the  reflections 
from  other  subsurface  objects  such  as  rocks  and  roots,  may  be  similar  to  those  from  a  mine. 

The  system  approach  taken  may  overcome  many  of  these  obstacles.  Initial  experiments 
were  performed  using  a  bistatic  forward  scattering  mode  as  described  in  QPR  No.  55  [1]. 
Since  the  oblique  forward  scattering  from  mine  geometries  is  very  high  compared  to 
backscattering,  a  much  larger  target  return  is  obtained  in  the  forward  direction.  More 
importantly,  the  soil  reflection  coefficient  decreases  dramatically  for  oblique  forward  scattering 
angles  near  the  Brewster  angle  as  discussed  in  [1].  Secondly,  the  use  of  a  high  resolution 
focused  array  will  allow  imaging  of  the  mine  itself  in  order  to  distinguish  it  from  false  targets. 
Initial  experiments  are  being  performed  at  L-band  at  a  frequency  of  1550  MHz  which  will  give 
resolution  on  the  order  of  the  size  of  a  mine.  High  resolution  would  be  obtained  at  S-band  (e.g. 
3500  MHz)  so  that  several  resolution  cells  would  exist  on  the  mine. 

This  section  will  report  on  initial  results  obtained  to  achieve  these  goals.  Included  are 
the  first  experimental  synthetic  aperture  array  image  of  a  buried  mine-like  target  and 
experimental  results  which  verify  the  forward  scattering  concept. 


2.2  Initial  Experimental  Results 


2.2.1  Data  Collection  and  Array  Processing  Software 

The  software  consists  of  two  primary  programs.  The  first  is  MINCOLLECT  which 
reads  and  stores  the  sampled  I  and  Q  channels  from  the  HP  54120T  sampling  scope  for  either 
single  event  experiments  or  for  each  of  the  16  element  channels  in  a  synthetic  aperture  array 
experiment.  This  program  completely  automates  the  experimental  data  collection  since  the  data 
is  cataloged  and  can  be  recalled  and  plotted  on  the  computer  CRT  screen  at  any  time.  The 
second  program  is  MINEPROC  which  is  the  synthetic  aperture  array  processing  program.  This 
program  gets  the  synthetic  aperture  array  data  for  16  elements  from  MINCOLLECT.  A  data 
range  bin  is  chosen  by  the  user  and  a  focusing  range  is  also  chosen.  The  program  then  digitally 
beamforms  in  angle  for  100  beam  positions  and  the  results  are  scaled  and  plotted  on  the  CRT  as 
response  vs.  angle.  Both  programs  have  been  tested  with  real  data  as  will  be  discussed  below. 

2.2.2  Propagation  Experiments 

The  purpose  of  these  experiments  is  to  determine  the  reflection  and  scattering 
characteristics  of  soil  and  targets  for  various  angles  of  incidence  for  both  backscattering  and 
forward  scattering  geometries.  The  experimental  setup  allows  measurement  of  backscattering 
for  angles  of  incidence  0°  to  90'  with  the  transmit  and  receive  antennas  mounted  adjacent  to  each 
other.  For  forward  scattering  measurements,  the  receive  antenna  is  oriented  to  the  angle  of  the 
forward  reflected  ray,  while  the  transmit  antenna  is  pointed  to  the  desired  angle  of  incidence. 
The  antenna  heights  and  separations  are  adjusted  to  insure  that  they  are  aimed  to  a  common 
midpoint  on  the  sand  surface. 

The  experimental  procedure  is  to  obtain  the  level  of  the  incident  wave  by  measuring  the 
direct  path  signal  with  the  antennas  pointed  at  each  other.  From  knowledge  of  this  level  and  of 
the  antenna  separation,  the  reflection  measurements  can  be  properly  normalized  to  compute 
reflection  coefficients.  A  normal  incidence  measurement  is  then  taken  at  a  known  height  and  the 
voltage  reflection  coefficient  is  computed.  From  this  measurement  the  dielectric  constant  is 
determined,  and  the  propagation  velocity  and  wavelength  in  the  sand  are  computed.  This 
measurement  is  necessary  since  the  moisture  content  of  the  sand  varies  on  a  daily  basis  causing 
variations  in  the  dielectric  constant  and  signal  attenuation. 

Four  types  of  propagation  measurements  have  been  taken  and  are  discussed  below.  The 
target  used  was  a  13"  x  17"  metal  plate.  The  nominal  pulse  width  is  8  nsec  (no  range  resolution 


is  needed  for  these  measurements).  Vertical  polarization  was  used  for  all  measurements. 
Horizontal  polarization  measurements  will  be  performed  later  for  comparison. 

a)  Interference  Measurements  (Vertical  Incidence! 

An  interesting  phenomenon  has  been  observed  when  the  target  is  buried  at  different 
depths  and  the  reflected  signal  is  measured.  The  received  signal  undergoes  a  cyclic  variation  in 
amplitude  due  to  interference  between  the  reflection  from  the  target  and  the  reflection  from  the 
sand,  as  described  in  Stratton  [9].  The  depth  of  maximum  signal  will  depend  upon  the  dielectric 
constant  of  the  medium.  These  data  are  now  being  analyzed  and  will  be  compared  with  an 
analytical  model.  This  interference  might  be  utilized  to  give  further  information  about  the  target 
such  as  depth  or  composition. 

Vertical  incidence  measurements  have  also  been  used  to  determine  dielectric  constant.  A 
typical  measured  value  of  the  voltage  reflection  coefficient  is  .37  which  gives  a  dielectric 
constant  of  4.73  and  a  Brewster  angle  (angle  of  incidence  for  zero  surface  reflection)  of  65.3°. 

b)  Backscattering  Measurements 

Preliminary  backscattering  measurements  have  been  done  with  the  metal  plate  target  at  a 
depth  of  2"  and  an  angle  of  incidence  of  45*.  The  measured  voltage  reflection  coefficient  is 
.108  compared  to  a  measured  forward  scattering  reflection  coefficient  approaching  unity. 
Although  there  was  no  measurable  backscattering  from  the  flat  sand  surface  alone, 
backscattering  from  a  rough  surface  could  easily  be  as  high  or  higher  than  that  from  the  target. 
Additional  measurements  will  be  performed  at  different  angles  to  further  characterize 
backscattering  properties. 

c)  Reflection  Coefficient  of  Sand  vs.  Angle  of  Incidence 

Because  of  the  utility  of  the  Brewster  angle  concept  for  vertical  polarization,  the  voltage 
reflection  coefficient  of  the  sand  was  measured  in  the  forward  scattering  setup  for  varying  angle 
of  incidence.  As  expected  its  value  decreased  from  about  .4  at  near  vertical  incidence  to  a 
minimum  of  .16  as  angle  increased  where  upon  it  again  began  to  increase  for  angles 
approaching  grazing.  The  minimum  occurs  near  an  angle  of  45*  rather  than  the  theoretical 
Brewster  angle,  probably  due  to  variations  in  the  angle  of  the  sand  surface  and  errors  in 
pointing  the  antennas.  Further  experiments  will  be  performed  with  these  errors  minimized. 


d) 


Forward  Scattering  from  Buried  Target 


The  previous  measurement  of  surface  reflections  for  varying  angle  of  incidence  was 
repeated  with  the  metal  plate  buried  at  a  depth  of  2”.  The  total  reflection  coefficient  (target  plus 
sand)  varies  from  near  unity  of  vertical  incidence  to  .59  at  60°  indicating  a  much  larger  forward 
scattering  amplitude  than  backscattering.  When  these  results  are  combined  with  the  measured 
sand  surface  reflection  coefficients,  a  signal-to-clutter  ratio  curve  can  be  computed  as  shown  in 
Figure  3.1.  Plotted  are  the  measured  values  of  the  forward  scattering  reflection  with  the  target 
present  normalized  by  those  with  the  target  absent.  The  peak  value  occurs  at  near  45°  with  a 
ratio  of  17.4  dB.  This  high  level  demonstrates  the  feasibility  of  detecting  buried  targets  in  a 
forward  scattering  mode. 

e)  Focused  Array  Imaging  Experiment 

The  first  experiment  to  detect  and  localize  a  buried  target  was  performed  using  a  16 
element  synthetic  aperture  focused  array.  The  transmit  antenna  and  receiving  array  were 
separated  by  125"  with  heights  of  36"  and  pointed  at  angles  of  incidence  and  reflection  of  60°. 
The  metal  plate  was  buried  at  a  depth  of  2"  at  the  midway  point  on  the  sand  between  the 
antennas.  The  single  receive  array  element  was  moved  to  16  sequential  positions  4"  apart  (1/2 
wavelength)  on  a  horizontal  line  perpendicular  to  the  plane  of  incidence.  For  each  element 
position,  the  I  and  Q  channel  signals  were  sampled  and  stored  by  the  sampling  scope  and  the 
data  were  read  by  the  program  MINCOLLECT  described  above.  The  synthetic  aperture  array 
processing  program  MINEPROC  was  then  used  to  focus  the  array  at  the  target  range  (1.86  m) 
and  to  scan  the  beam  to  100  sequential  positions  in  horizontal  angle  at  the  target  range  and 
compute  the  response.  The  resulting  image  of  the  metal  target  is  shown  in  Figure  3.2.  The 
peak  response  is  at  the  target  location  and  has  a  width  approximately  equal  to  the  angular  width 
of  the  target  The  sidelobe  response  of  the  array  is  seen  to  be  very  good,  on  the  order  of  about  - 
13  dB,  the  theoretical  level  of  a  uniform  array  (no  amplitude  tapering  was  used  for  this  initial 
image  although  it  can  easily  be  incorporated  into  the  program).  Sequential  images  with  the  array 
focused  at  ranges  other  than  the  target  range  show  a  gradual  decrease  in  peak  amplitude. 

The  relative  peak  voltage  amplitude  of  the  target  response  is  1261.  This  is  to  be 
compared  to  the  peak  response  of  the  image  when  the  target  is  absent,  i.e.,  of  the  sand  alone. 
This  image  is  shown  in  Figure  3.3  where  the  peak  amplitude  has  decreased  to  281.  Thus  the 
peak  response  with  the  target  present  is  higher  than  that  of  the  sand  alone  by  13.8  dB. 
Therefore  the  peak  response  to  the  sand  alone  is  less  than  the  sidelobe  response  to  the  target. 


These  preliminary  results  indicate  a  favorable  chance  for  success  in  imaging  of 
subsurface  mines  with  a  focused  synthetic  aperture  array.  Additional  experimental  plans  include 
the  use  of  target  of  varying  dielectric  constant  including  real  mine  casings,  experiments  in  a 
moist  sand  medium,  further  forward  and  backscattering  experiments  for  additional  angles  and 
ranges,  and  continuation  of  array  imaging  experiments  under  these  conditions.  The  use  of  a  1 
nsec  pulse  and  the  use  of  S-band  hardware  can  allow  even  higher  resolution  imaging. 


2.3  Mine  Detection  System 


The  system  concept  proposed  is  the  use  of  a  bistatic  radar  system  operating  in  the 
frequency  range  of  L-band  to  S-band.  The  lower  frequencies  are  more  useful  for  minimizing 
attenuation,  while  the  higher  frequency  range  yields  more  desirable  resolution  for  imaging  the 
mine  itself.  After  a  review  of  the  system  requirements  it  was  decided  to  conduct  first  year 
experiments  with  a  4  nsec  pulse  at  L-band  initially,  despite  inadequate  range  resolution  for  a 
shallow  subsurface  target.  The  receiver  consists  of  a  horizontal  line  array  of  equally  spaced 
elements  with  the  image  space  located  in  the  Fresnel  region  of  its  field  of  view.  The 
transmitting  antenna  illuminates  the  ground  at  an  oblique  angle,  and  the  receiving  array 
measures  the  forward  scatter  from  the  medium.  Each  array  element  then  samples  and  stores  the 
time  sequence  from  each  pulse.  Digital  beamforming  is  performed  in  the  system  computer  by 
focusing  and  scanning  the  array  in  two  dimensions  (transverse  and  depth).  Imaging  in  the  third 
dimension  (longitudinal)  may  be  assisted  by  forward  motion  of  the  supporting  vehicle. 

2.3.1  Basic  System  Features 

The  initial  experiment  is  performed  at  1.5  Ghz  using  a  horn  antenna  at  the  transmitter, 
and  dipole/comer-reflector  combination  antenna  at  the  single  receiver.  A  500  Mhz  scope  is  used 
with  the  receiver  to  view  the  envelope  of  the  received  waveform  provided  by  an  envelope 
detector.  A  pulse  generator  is  used  as  the  trigger  for  a  diode  switch  that  generates  a  1  nsec  pulse 
modulation  of  the  source  frequency.  The  experiment  is  performed  in  a  1-foot  deep  sandbox  of 
dimensions  8  ft  x  16  ft.  The  experimental  setup  is  based  on  a  20  Ghz  bandwidth  HP  54120-T 
sampling  scope.  This  allows  sampling  and  storage  of  the  entire  broadband  waveform.  The 
sampling  scope  interfaces  with  the  HP  9816S  computer  which  can  input,  store  and  display  the 
waveform  data  and  perform  necessary  computations.  The  data  can  be  transferred  to  the  Moore 
School  laboratory  for  further  processing  on  the  compatible  HP  310  or  on  a  VAX  system. 

2.3.2  A  Method  of  Eliminating  Ground  Surface  Reflections 

One  of  the  biggest  acknowledged  problems  in  the  use  of  subsurface  radar  for  mine 
detection  is  the  high  amplitude  reflection  from  the  air-ground  interface.  For  a  normally  incident 
plane  wave,  the  ground-reflected  amplitude  may  be  as  high  as  50%  of  that  of  the  incident  wave, 
while  the  reflection  from  a  mine-like  target  may  be  much  less.  This  would  make  detection  of  the 
target  impossible  if  it  were  within  the  same  pulse  return  as  the  surface  interface.  This  makes 
necessary  the  use  of  a  very  short  pulse  for  a  shallow  target  in  order  to  distinguish  target 
response  from  the  surface  interface  reflection.  This  imposes  severe  requirements  on  the  radar 


system  with  the  conflicting  need  for  a  wide  bandwidth,  but  low  frequency  to  minimize 
attenuation.  Even  if  this  is  done,  the  signal  loss  from  surface  reflection  may  not  be  affordable. 

A  potential  solution  to  this  problem  is  proposed;  the  solution  utilizes  the  polarizing  angle 
or  Brewster  angle  to  achieve  100%  transmission  of  incident  wave  into  the  earth.  If  the  incident 
wave  is  polarized  with  its  electric  field  parallel  to  the  plane  of  incidence  (corresponding  to  a 
vertically  polarized  antenna),  and  if  the  wave  is  incident  at  the  air-ground  Brewster  angle,  then 
there  will  be  no  surface  reflections.  The  wave  transmitted  into  the  ground  will,  upon  reflection 
from  the  target  (assuming  geometical  optics  reflection),  be  incident  on  the  ground-air  interface  at 
the  ground-air  Brewster  angle,  and  thus  be  completely  transmitted  into  the  air.  Thus,  the 
reflection  at  both  interfaces,  air-ground  and  ground-air,  are  eliminated  with  100%  transmission 
back  to  the  receiver. 

A  possible  scheme  for  realization  of  this  type  of  propagation  geometry  is  shown  in 
Figure  3.4.  For  dry  soil  with  a  dielectric  constant  of  3,  the  Brewster  angle  is  60”  for  the 
ground-air  interface.  The  transmitting  antenna  should  be  rotatable  and  translatable  to  achieve 
this  and  to  adjust  for  varying  Brewster  angles. 

There  are  many  possible  variations  of  geometry  for  the  receiving  array.  The  most 
obvious  one  is  with  the  array  perpendicular  to  the  plane  of  incidence.  Both  transmit  and  receive 
antennas  must  be  vertically  polarized.  The  receiving  array  should  also  be  capable  of  rotation 
and  translation  in  order  to  optimize  the  return  from  various  portions  of  the  subsurface  area  to  be 
imaged. 

Another  possible  geometry  of  the  receiving  array  is  in  the  plane  of  incidence.  In  this 
case,  different  array  elements  would  have  varying  responses  to  reflections  from  targets  in 
different  regions. 

Another  advantage  of  this  type  of  propagation  geometry  is  that  the  receiving  array  would 
be  located  in  the  region  of  forward  scattering  of  the  target,  a  region  of  very  high  intensity.  An 
oblique  forward-looking  monostatic  antenna  cannot  be  relied  upon  because  of  the  low 
backscattering  cross  sections  of  typical  targets  at  oblique  angles. 

There  are  many  other  possibilities  as  well  as  problems  using  this  concept  for  mine 
detection.  A  major  consideration  is  the  critical  angle  at  which  total  internal  reflection  of  the 
target  energy  occurs  at  the  ground-air  interface.  This  angle  is  always  greater  than  the  Brewster 
angle,  but  the  angles  become  closer  as  the  dielectric  constant  of  the  ground  increases. 

Another  factor  is  the  finite  conductivity  of  the  ground  when  there  is  some  moisture  content.  The 
wavenumber  then  becomes  complex,  and  reflection  and  refraction  depend  upon  both  dielectric 
constant  and  the  conductivity  of  the  medium. 

The  measurements  which  should  be  performed  in  this  experiment  include: 


Figure  3.4.  Brewster  Angle  Propagation  Geometry 


( 1 )  reflection  coefficient  measurements  from  air-ground  interface  for  varying  angle  of 
incidence,  and  various  soil  types  and  moisture  contents; 

(2)  free-space  scattering  measurements  for  various  target  dielectrics,  sizes  and  shapes 
for  varying  angle  of  incidence; 

(3)  scattering  properties  of  the  targets  of  (2)  when  placed  in  the  medium  of  (1); 

The  above  measurements  should  be  performed  for  varying  depths  and  horizontal 
displacements  of  the  targets,  as  well  as  varying  polarization,  transmit  antenna  directivity, 
and  height  of  T/R  antennas. 

2.3.3  Experimental  Results 

a)  System  Approach 

The  goal  of  this  program  is  to  study  the  application  of  a  short  pulse,  high  angular 
resolution  array  system,  focused  in  the  near  field,  to  the  detection  and  imaging  of  shallow 
subsurface  land  mines.  The  method  of  implementation  is  to  be  synthetic  aperture  digital 
beamforming,  in  which  the  array  element  outputs  are  sequentially  sampled,  and  the 
beamforming,  signal  processing  and  imaging  are  performed  in  the  system  computer. 

The  major  obstacles  to  achieving  these  goals  are  the  varying  dielectric  constant  and 
conductivity  of  typical  soil  types.  Land  mines  are  typically  only  a  few  centimeters  below  the 
soil  surface.  The  relatively  high  dielectric  constant  of  the  soil  will  cause  a  large  ground  surface 
reflection  which  may  be  larger  than  the  return  from  a  target,  and  both  returns  may  be  in  the  same 
range  bin;  even  a  pulse  as  short  as  one  nanosecond  has  a  range  resolution  in  sand  of  10  cm  for 
backscatter,  20  cm  for  forward  scatter,  and  probably  could  not  resolve  a  shallow  target  from  the 
surface  reflection.  The  finite  conductivity  of  the  soil  may  cause  a  large  attenuation  of  the  signal 
for  high  frequencies,  as  well  as  phase  variations  in  the  signal  return.  Also,  target  dielectric 
constant  may  be  similar  to  that  of  the  soil  making  detection  difficult.  Finally,  the  reflections 
from  other  subsurface  objects  such  as  rocks  and  roots,  may  be  similar  to  those  from  a  mine. 

The  system  approach  taken  overcomes  many  of  these  obstacles.  Initial  experiments 
were  performed  using  a  bistatic  forward  scattering  mode  as  described  in  Valley  Forge  Research 
Center  Quarterly  Progress  Report  No.  55  [1].  Since  the  oblique  forward  scattering  from  mine 
geometries  is  very  high  compared  to  backscattering,  a  much  larger  target  return  was  obtained  in 
forward  direction.  More  importantly,  the  soil  reflection  coefficient  decreases  dramatically  for 
oblique  forward  scattering  angles  near  the  Brewster  angle  as  discussed  in  [1].  Secondly,  the 
use  of  a  high  resolution  focused  array  will  allow  imaging  of  the  mine  itself  in  order  to 
distinguish  it  from  false  targets.  Initial  experiments  were  performed  at  L-band  at  a  frequency  of 


1550  MHz,  which  gives  resolution  on  the  order  of  the  size  of  a  mine.  Higher  resolution  would 
be  obtained  at  S-band  (e.g.,  3.5  GHz)  so  that  several  resolution  cells  would  exist  on  the  mine. 

This  section  reports  on  the  initial  results  obtained  to  achieve  these  goals.  Included  are 
the  first  experimental  synthetic  aperture  array  image  of  a  buried  mine-like  target  and 
experimental  results  which  verify  the  forward  scattering  concept. 

b)  Data  Collection  and  Array  Processing  Software 

The  software  consists  of  two  primary  programs.  The  first  is  MINCOLLECT  which 
reads  and  stores  the  sampled  I  and  Q  channels  from  the  HP  54120T  sampling  scope  for  either 
single  event  experiments  or  for  each  of  the  16-element  channels  in  a  synthetic  aperture  array 
experiment.  This  program  completely  automates  the  experimental  data  collection  since  the  data 
is  catalogued  and  can  be  recalled  and  plotted  on  the  computer  CRT  screen  at  any  time.  The 
second  program  is  MINEPROC  which  is  the  synthetic  aperture  array  processing  program.  This 
program  gets  the  synthetic  aperture  array  data  for  16  elements  form  MINCOLLECT.  A  data 
range  bin  is  chosen  by  the  user  and  a  focusing  range  is  also  chosen.  The  program  then  digitally 
beamforms  in  angle  for  100  beam  positions  and  the  results  are  scaled  and  plotted  on  the  CRT  as 
response  vs.  angle.  Both  programs  have  been  tested  with  the  real  data  as  will  be  discussed 
below. 

c)  Propagation  Experiments 

The  purpose  of  these  experiments  has  been  to  determine  the  reflection  and  scattering 
characteristics  of  soil  and  targets  for  various  angles  of  incidence  for  both  backscattering  and 
forward  scattering  geometries.  The  experimental  setup  allows  measurement  of  backscattering 
for  angles  of  incidence  0°  to  90’  with  the  transmit  and  receive  antennas  mounted  adjacent  to  each 
other.  For  forward  scattering  measurements,  the  receive  antenna  is  oriented  to  the  angle  of  the 
forward  reflected  ray,  while  the  transmit  antenna  is  pointed  to  the  desired  angle  of  incidence. 
The  antenna  heights  and  separations  are  adjusted  to  insure  that  they  are  aimed  to  a  common 
midpoint  on  the  sand  surface. 

The  experimental  procedure  is  to  obtain  the  level  of  the  incident  wave  by  measuring  the 
direct  path  signal  with  the  antennas  pointed  at  each  other.  From  knowledge  of  this  level  and  of 
the  antenna  separation,  the  reflection  measurements  can  be  properly  normalized  to  compute 
reflection  coefficients.  A  normal  incidence  measurement  is  then  taken  at  a  known  height  and  the 
voltage  reflection  coefficient  is  computed.  From  this  measurement  the  dielectric  constant  is 
determined,  and  the  propagation  velocity  and  wavelength  in  the  sand  are  computed.  This 
measurement  is  necessary  since  the  moisture  content  of  the  sand  varies  on  a  daily  basis  causing 
variations  in  the  dielectric  constant  and  signal  attenuation. 


Four  types  of  propagation  measurements  have  been  taken  and  are  discussed  below.  The 
target  used  was  a  13"  x  17"  metal  plate.  The  nominal  pulse  width  is  8  nsec  (no  range  resolution 
is  needed  for  these  measurements).  Vertical  polarization  was  used  for  all  measurements. 
Horizontal  polarization  measurements  will  be  performed  later  for  comparison. 

d)  Interference  Measurements  (Vertical  Incidence) 

An  interesting  phenomenon  has  been  observed  when  the  target  is  buried  at  different 
depths  and  the  reflected  signal  is  measured.  The  received  signal  undergoes  a  cyclic  variation  in 
amplitude  due  to  interference  between  the  reflection  from  the  target  and  the  reflection  from  the 
sand,  as  described  in  Stratton  [9].  The  depth  of  maximum  signal  will  depend  upon  the  dielectric 
constant  of  the  medium. 

Vertical  incidence  measurements  have  also  been  used  to  determine  dielectric  constant.  A 
typical  measured  value  of  the  voltage  reflection  coefficient  is  .37  which  gives  a  dielectric 
constant  of  4.73  and  a  Brewster  angle  (angle  of  incidence  for  zero  surface  reflection)  of  65.3°. 

e)  Backscattering  Measurements 

Preliminary  backscattering  measurements  have  been  done  with  the  metal  plate  target  at  a 
depth  of  2"  and  an  angle  of  incidence  of  45*.  The  measured  voltage  reflection  coefficient  is 
.108  compared  to  a  measured  forward  scattering  reflection  coefficient  approaching  unity. 
Although  there  was  no  measurable  backscattering  from  the  flat  sand  surface  alone, 
backscattering  from  a  rough  surface  could  easily  be  as  high  or  higher  than  that  from  the  target. 

f)  Reflection  Coefficient  of  Sand  vs.  Angle  of  Incidence 

Because  of  the  utility  of  the  Brewster  angle  concept  for  vertical  polarization,  the  voltage 
reflection  coefficient  of  the  sand  was  measured  in  the  forward  scattering  setup  for  varying  angle 
of  incidence.  As  expected  its  value  decreased  from  about  0.4  at  near  vertical  incidence  to  a 
minimum  of  0.16  as  angle  increased  whereupon  it  again  began  to  increase  for  angles 
approaching  grazing.  The  minimum  occurs  near  an  angle  of  45°  rather  than  the  theoretical 
Brewster  angle,  probably  due  to  variations  in  the  angle  of  the  sand  surface  and  errors  in 
pointing  the  antennas. 

g)  Forward  Scattering  from  Buried  Target 

The  previous  measurement  of  surface  reflections  for  varying  angle  of  incidence  was 
repeated  with  the  metal  plate  buried  at  a  depth  of  2".  The  total  reflection  coefficient  (target  plus 
sand)  varies  from  near  unity  at  vertical  incidence  to  0.59  at  60°  indicating  a  much  larger  forward 
scattering  amplitude  than  backscattering.  When  these  results  are  combined  with  the  measured 


sand  surface  reflection  coefficients,  a  signal-to-clutter  ratio  curve  can  be  computed  as  shown  in 
Figure  3.5.  Plotted  are  the  measured  values  of  the  forward  scattering  reflection  with  the  target 
present  normalized  by  those  with  the  target  absent.  The  peak  value  occurs  at  near  45°  with  a 
ratio  of  17.4  dB.  This  high  level  demonstrates  the  feasibility  of  detecting  buried  targets  in  a 
forward  scattering  mode. 

h)  Focused  Array  Imaging  Experiments 

The  first  experiment  to  detect  and  localize  a  buried  target  was  performed  using  a  16- 
element  synthetic  aperture  focused  array.  The  transmit  antenna  and  receiving  array  were 
separated  by  125"  with  heights  of  36"  and  pointed  at  angles  of  incidence  and  reflection  of  60°. 
The  metal  plate  was  buried  at  a  depth  of  2"  at  the  midway  point  on  the  sand  between  the 
antennas.  The  single  receive  array  element  was  moved  to  16  sequential  position  4"  apart  (1/2 
wavelength)  on  a  horizontal  line  perpendicular  to  the  plane  of  incidence.  For  each  element 
positions,  the  I  and  Q  channel  signals  were  sampled  and  stored  by  the  sampling  scope  and  the 
data  were  read  by  the  program  MINCOLLECT  described  above.  The  synthetic  aperture  array 
processing  program  MINEPROC  was  then  used  to  focus  the  array  at  the  target  range  (1.86  m) 
and  to  scan  the  beam  to  100  sequential  positions  in  horizontal  angle  at  the  target  range  and 
compute  the  response.  The  resulting  image  of  the  metal  target  is  shown  in  Figure  3.6.  The 
peak  response  is  at  the  target  location  and  has  a  width  approximately  equal  to  the  angular  width 
of  the  target.  The  sidelobe  response  of  the  array  is  seen  to  be  very  good,  on  the  order  of  about 
-13  dB,  the  theoretical  level  of  a  uniform  array  (no  amplitude  tapering  was  used  for  this  initial 
image  although  it  can  easily  be  incorporated  into  the  program).  Sequential  images  with  the  array 
focused  at  ranges  other  than  the  target  range  show  a  gradual  decrease  in  peak  amplitude. 

The  relative  peak  voltage  amplitude  of  the  target  response  is  1261.  This  is  to  be 
compared  to  the  peak  response  of  the  image  when  the  target  is  absent,  i.e.,  of  the  sand  alone. 
This  image  is  shown  in  Figure  3.7  where  the  peak  amplitude  has  decreased  to  281.  Thus  the 
peak  response  with  the  target  present  is  higher  than  that  of  the  sand  alone  by  13.8  dB. 
Therefore  the  peak  response  to  the  sand  alone  is  less  than  the  sidelobe  response  to  the  target. 

These  preliminary  results  indicate  a  favorable  chance  for  success  in  imaging  of 
subsurface  mines  with  a  focused  synthetic  aperture  array.  Future  experimental  plans  include  the 
use  of  targets  of  varying  dielectric  constant  including  real  mine  casings,  experiments  in  a  moist 
sand  medium,  further  forward  and  backscattering  experiments  for  additional  angles  and  ranges, 
and  additional  array  imaging  experiments  under  these  conditions.  A  1  nsec  pulse  capability 
gives  improved  range  resolution.  The  addition  of  S-band  capability  to  allow  higher  resolution 
imaging  is  suggested. 


Figure  3.5.  Ratio  of  forward  reflection  of  buried  metal  target  (T)  to 
sand  (S)  vs.  angle  of  incidence. 
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Figure. 3.6.  Buried  Metal  Plate. 
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Figure  3.7.  Sand  Alone. 


2.4  L-Band  Mine  Detection  Experiments 

The  outdoor  L-band  mine  detection  measurements  continued  after  the  completion  of  the 
indoor  mine  radar  cross-section  measurements.  The  objective  of  the  measurements  was  to 
verify  the  Brewster  angle  effect  on  reduction  of  reflected  energy  in  a  bistatic  oblique  forward 
scattering  mode.  It  had  been  found  during  the  indoor  RCS  measurements  that  the  direct  path 
radiation  between  antennas  was  large  enough  to  cause  error  in  the  measurements.  A  centrally 
located  absorber  was  used  between  the  transmit  and  receive  antennas  to  successfully  eliminate 
the  direct  path  interference.  This  same  approach  was  used  in  the  outdoor  measurements.  A 
microwave  absorber  barrier  was  suspended  between  the  two  antennas  to  block  the  direct  path 
radiation  but  allow  the  ground  reflected  energy  to  propagate  under  the  barrier. 

Reflection  coefficient  measurements  were  made  on  the  sand  by  varying  the  angle  of 
incidence  between  40°  and  75°.  A  distinct  Brewster  angle  effect  was  seen  in  all  the 
measurements,  and  the  reflection  coefficient  showed  a  noticeable  decrease  over  the  range  of  45° 
to  60°  angle  of  incidence.  The  minimum  measured  value  of  the  reflection  coefficient  was  0.068 
at  about  60°  compared  to  a  value  of  0.231  at  45°. 

Another  method  was  discovered  to  reduce  surface  reflections  in  bistatic  forward 
scattering.  Since  the  ground  surface  was  in  proximity  to  the  transmit  and  receive  antennas,  the 
reflected  energy  distribution  was  approximately  the  geometrical  projection  of  the  transmit  horn 
aperture.  Since  the  measured  scattering  from  the  mines  had  been  found  to  be  wide  angle,  it  was 
actually  possible  to  place  the  receive  antenna  such  that  it  was  just  outside  of  the  surface-reflected 
tube  of  energy  from  the  transmitter  but  still  near  the  maximum  of  the  mine  scattering  pattern. 

This  technique  was  used  to  make  synthetic  aperture  measurements  of  the  bistatic 
forward  scattered  power  from  buried  mines.  The  angle  of  incidence  used  was  45°.  The  focal 
distance  of  the  array  was  1.8  m.  The  results  are  illustrated  in  Figures  3.8  -  3.10.  Figure  3.8 
shows  the  image  as  a  function  of  azimuth  angle  of  the  surface  reflection  with  no  buried  mine. 
The  peak  level  is  -14  dB  as  shown.  The  square  plastic  M-19  anti-tank  mine  was  then  buried  at  a 
depth  of  7.5  cm  and  the  image  of  Figure  3.9  was  obtained  showing  a  peak  of  -9  dB,  an  increase 
of  5  dB.  This  is  a  large  increase  in  power  for  this  type  of  mine  since  its  dielectric  constant  is 
close  to  that  of  sand  and  it  consequently  has  a  low  reflectivity.  The  same  procedure  was  used 
for  the  metal  anti-tank  mine  buried  at  a  depth  of  5  cm.  In  this  case  (see  Figure  3.10),  the  peak 
power  is  -2.5  dB,  an  increase  of  11.5  dB  over  the  sand  surface  reflection  alone.  These 
experiments  have  demonstrated  the  ability  of  our  technique  to  detect  both  types  of  mine  with  a 
high  signal-to-clutter  ratio. 
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Figure  3.8.  SAR  Forward  Scatter.  Sand  Alone. 

AOI  =  45*.  Focal  Distance  =  1.8  m. 


Figure  3.9.  SAR  Forward  Scatter.  Plastic  Mine. 

Depth  =  7.5  cm.  AOI  =  45'.  Focal  Distance  =  1.8  m. 
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Figure  3.10 


SAR  Forward  Scatter.  Metal  Mine. 

Depth  =  5  cm.  AOI  =  45’.  Focal  Distance  =  1.8  m. 


2 . 5  S-Band  Experiments 


Experiments  were  at  S-band  using  a  frequency  of  3.5  GHz.  The  experimental  hardware 
is  described  later  in  this  article.  The  purpose  of  increasing  the  frequency  to  S-band  is  to 
improve  the  resolution  to  allow  imaging  of  the  mine  for  identification  purposes. 

Initial  measurements  have  been  made  of  the  forward  reflection  from  the  sand  surface; 
these  measurements  have  also  shown  a  significant  decrease  in  reflection  at  the  Brewster  angle 
(55*  -  60°).  Another  valuable  feature  of  the  S-band  setup  is  the  narrow  pulse  width  of  1  nsec 
achievable  with  the  diode  switch,  compared  to  about  10  nsec  at  L-band.  The  greatest  advantage 
of  this  is  the  ability  to  discriminate  against  the  direct  path  radiation  because  of  its  path  length 
difference  from  the  reflected  path. 

Measurements  have  also  been  made  of  the  reflection  from  buried  mines,  both  with  a 
single  receiver,  and  with  the  synthetic  aperture  array.  The  1.5  m  size  of  the  array  aperture  is  the 
same  as  at  L-band  with  the  same  number  of  elements,  16,  now  spaced  at  about  1  wavelength. 
Since  the  field  of  view  of  the  array  is  only  about  45°,  the  grating  lobe  at  90°  will  not  be  visible. 
These  data  have  been  analyzed  and  are  discussed  in  following  sections. 

2.6  Mine  Scattering 

The  purpose  of  this  task  is  to  determine  the  forward  reflection  characteristics  of  anti-tank 
mines  by  analysis  of  experimental  data.  Experiments  have  been  performed  with  metal  and 
plastic  mines  in  both  free-space  and  in  soil.  A  horizontal  line  array  was  also  used  with  buried 
mines  to  determine  their  detectability  in  an  oblique  forward  scattering  mode  near  the  Brewster 
angle.  The  measurements  were  made  at  L-band  (1.55  Ghz)  and  S-band  (3.5  Ghz).  The  soil 
measurements  used  sand  with  varying  amounts  of  moisture. 

2.6.1  Radar  Cross-Section  Measurements 

Measurements  of  the  radar  cross-section  as  a  function  of  angle  of  incidence  were 
performed  at  L-band  for  both  metal  and  plastic  anti-tank  mines.  The  purpose  was  to  determine 
if  the  mines  had  a  significant  forward  reflected  signal  when  the  transmitter  was  located  at  an 
oblique  angle  of  incidence  and  the  receiver  was  at  the  opposite  angle  of  reflection  with  the  mine 
in  between,  as  in  the  geometry  of  Figure  3.11.  The  measurements  were  made  at  vertical 
polarization  (electric  field  in  the  plane  of  incidence)  which  is  the  system  polarization  to  be  used 
for  the  Brewster  angle  effect.  Monostatic  backscattering  measurements  were  also  made  for 


Figure  3.1 1(a).  Support  Mounting  of  Mine,  for  Experiment. 
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(a)  Forward  scattering  geometry. 


<b>  Backscattering  geometry. 


Figure  3.11(b).  Scattering  Geometry. 


comparison.  The  measurements  were  made  in  free  space  using  a  narrow  pulse  so  that  all  other 
reflections  were  eliminated. 

The  results  show  an  almost  constant  return  from  the  mines  as  a  function  of  angle  of 
incidence  for  the  forward  scattering  mode  out  to  wide  angles.  Since  the  angle  of  incidence  on 
the  mines  would  be  smaller  when  they  are  buried  due  to  refraction,  the  results  show  that  the 
oblique  forward  reflection  from  a  mine  would  be  the  same  as  normal  incidence.  Typical 
forward  scattering  results  for  the  metal  M15  mine  and  the  plastic  M19  mine  are  shown  in 
Figures  3.12  and  3.13.  Figure  3.12  is  the  RCS  at  an  L-band  frequency  of  1550  MHz  as  a 
function  of  angle  of  incidence  on  the  plane  of  the  mine  (same  as  angle  of  reflection)  and  is 
almost  constant  at  its  maximum  out  to  30‘ .  Since  the  angle  of  refraction  in  soil  is  less  than  30° 
for  angles  of  incidence  of  less  than  60‘ ,  the  forward  reflection  from  the  mine  would  be  almost 
constant  at  all  oblique  angles  of  incidence.  Figure  3.13  is  the  corresponding  result  for  the  Ml 5 
and  shown  similar  behavior. 

These  are  in  contrast  to  backscattering  shown  in  Figures  3.14  and  3.15  in  which  the 
transmitter  and  receiver  are  at  the  same  location.  The  reflected  backscatter  signal  decreases 
rapidly  as  the  angle  of  incidence  on  the  plane  of  the  mine  increases  beyond  5°.  Similar  results 
were  obtained  at  a  frequency  of  f=1300  MHz  with  both  mines. 

The  S-band  forward  scattering  results  at  3500  MHz  are  given  in  Figures  3.16  and  3.17. 
Figure  3.16  shows  the  M15  scattering  decreasing  slightly  beyond  normal  incidence,  and  then 
continually  increasing  to  its  maximum  level  at  wide  angles.  The  Ml 9  result  in  Figure  3.17 
shows  a  relatively  constant  level  out  to  wide  angles.  The  backscattering  results  of  Figures  3.18 
and  3.19  demonstrate  a  series  of  rapid  variations  between  maxima  and  minima.  The  case  of  the 
M19  shown  a  main  backscattering  lobe  only  5°  wide. 

These  results  demonstrate  the  realiability  and  consistency  of  high  scattering  levels  at 
wide  angles  in  a  bistatic  forward  scattering  mode.  Measurements  at  horizontal  polarization 
showed  similarly  high  wide  angle  scattering. 

2.6.2  Buried  Mine  Measurements 

The  mine  forward  scattering  measurements  were  also  performed  in  dry  and  moist  sand 
for  fixed  angle  of  incidence  and  varying  mine  depth.  A  summary  of  some  of  the  results  is  given 
in  Table  1.  The  ratio  of  the  signal -to-clutter  of  the  mine  reflection  to  the  soil  reflection  is  given 
as  S/C.  The  M15  metal  mine  has  an  expected  high  forward  scattering  level  at  both  frequencies. 
The  plastic  M19  has  a  lower  return  but  would  still  be  detectable,  especially  when  an  array  is 
used. 
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Figure  3.12.  Forward  Scattering  RCS  of  Metal  Mine  at  L-Band. 


SCATTERING  ANGLE  (deg) 

FILE**Di  R289020B.  RCS 
FORWARD  SCATTERING 

M19  PLASTIC  MINE  f*1550  MHz  E-POL  R=2.24  m 


Figure  3.13.  Forward  Scattering  RCS  of  Plastic  Mine  at  L-Band. 
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Figure  3.14.  Backscattering  RCS  of  Metal  Mine  at  L-Band. 
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Figure  3.15.  Backscattering  RCS  of  Plastic  Mine  at  L-Band. 
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Figure  3.16.  Forward  Scattering  RCS  of  Metal  Mine  at  S-Band. 


SCATTERING  ANGLE  (deg) 
FILE-Dj  R29001 17. RC6 
FORWARD  SCAT  I* EE. Ft  1  fs| C3 


M19  f-3500  MHz  E-POL  R=2.44  m 


Figure  3.17.  Forward  Scattering  RCS  of  Plastic  Mine  at  S-Band 
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Figure  3.19.  Backscattering  RCS  of  Plastic  Mine  at  S-Band. 
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Table  1.  Peak  Reflection  Levels 


L-Band 

(1.55  GHz)  0%  water 

Mine 

R 

AOI 

d 

S/C 

M15 

1.8m 

45°  5  cm 

12  dB 

M19 

1.8m 

45°  7.5  cm 

5dB 

S-Band 

(3.5  GHz) 

0%  water 

Mine 

R 

AOI 

d 

S/C 

M15 

1.2m 

55°  3.8  cm 

20  dB 

M19 

1.2m 

55°  3.5  cm 

5  dB 

S-Band 

(3.5  GHz) 

10%  water 

Mine 

R 

AOI 

d 

S/C 

M19 

1.2m 

55'  3.8  cm 

4  dB 

2.6.3  Focused  Array  Measurements 

Focused  array  measurements  were  performed  using  the  forward  scattering  horizontal 
line  array  setup.  The  arrays  were  synthesized  by  sequentially  moving  a  single  antenna  to  each 
of  the  desired  array  element  locations  and  sampling  the  element  I  and  Q  outputs.  The  array  was 
then  digitally  focused  and  scanned  to  form  the  images.  Experiments  were  performed  using  a  16 
element  L-band  array  and  a  22  element  S-band  array.  The  array  response  results  show 
detectibility  of  both  metal  and  plastic  buried  mines  at  L-band  but  no  shape  information  due  to  the 
long  wavelength  (20  cm).  In  contrast,  the  S-band  array  data  shown  a  distince  variation  in  level 
while  scanning  over  the  surfaces  of  both  metal  M15  and  plastic  M19  buried  mines.  This  is  an 
initial  demonstration  of  the  ability  to  obtain  shape  information  for  buried  mines  at  S-band 
frequencies. 

These  S-band  measurements  of  anti-tank  mines  buried  in  sand  were  made  using  a 
bistatic  S-band  radar  system  having  a  wide  beam  standard  gain  horn  as  the  illuminating  antenna 
and  a  22  element  horizontal  line  array  as  the  receiving  antenna.  The  system  used  the  oblique 


forward  scattering  geometry  with  the  transmit  and  receive  antennas  at  the  55°  Brewster  angles  of 
incidence  and  reflection,  respectively.  The  antennas  were  at  a  height  above  the  sand  surface  of 
.7  m  with  a  separation  of  2.0  m.  The  targets  were  individually  located  centrally  between  the 
transmit  and  receive  antennas,  and  the  baseband  I  and  Q  element  data  was  sampled  by  the  20 
GHz  HP54120T  sampling  oscilloscope.  Focused  array  digital  beamforming  was  then 
performed  in  the  cross-range  dimension  with  focusing  depth  as  a  parameter.  The  array  was 
scanned  to  100  points  in  cross-range  (transverse  horizontal)  over  aim  distance,  at  each  of  25 
focal  depths  in  4  cm  increments  varying  from  .2  m  above  the  surface  of  the  ground  to  .8  m 
below  the  surface.  The  resulting  data  was  then  converted  to  power,  quantized  into  16  levels, 
and  a  half-tone  gray-level  image  was  generated  for  each  target  with  power  decreasing  from 
white  to  black. 

The  image  in  Figure  3.20  shows  the  response  of  the  array  to  the  sand  in  the  absence 
of  any  target,  where  the  image  is  normalized  to  its  own  peak  (voltage)  amplitude  of  23  (in 
relative  units).  The  image  in  Figure  3.21  is  that  of  a  1  in  wide  metal  ruler  oriented 
longitudinally  on  the  surface  of  the  sand  in  order  to  test  the  resolution  capability  of  the 
array.  The  data  is  normalized  to  its  peak  amplitude  value  of  133,  thus  the  sand  surface 
reflection  is  below  the  lowest  quantization  level  and  is  invisible.  From  the  actual  amplitude 
data  the  resolution  is  ascertained  to  be  7.5  cm.  The  theoretical  depth  of  field,  (depth 
resolution)  appears  to  be  near  its  theoretical  value  of  30  cm.  The  peak  intensity  of  the 
image  also  appears  correctly  at  the  location  of  the  surface.  The  image  of  Figure  3.22  is  that 
of  an  M15  metal  mine  buried  at  a  depth  of  7.5  cm  and  again  normalized  to  its  peak  of  139, 
with  the  surface  reflection  also  below  the  lowest  quantization  level  and  appearing  black.  In 
this  image,  there  is  notable  detail  across  the  mine  horizontal  dimension  with  the  largest 
reflection  for  the  raised  central  pressure  plate  and  two  secondary  reflection  lobes  from  the 
annulus  surrounding  the  pressure  plate.  The  highest  image  intensity  also  occurs  at  the 
correct  depth  of  7.5  cm  and  occurs  significantly  lower  in  position  than  the  peak  of  the  ruler 
image.  Some  smearing  of  the  image  occurs  above  the  surface  since  the  depth  of  field  is 
large  enough  to  capture  some  of  the  very  high  reflections  from  the  mine  surface  even  while 
it  is  focused  above  the  ground.  This  phenomenon  will  not  occur  with  weaker  targets,  such 
as  the  M19  plastic  mine  also  at  a  depth  of  7.5  cm,  shown  in  Figure  3.23.  The  peak  relative 
amplitude  of  the  mine  is  47,  and  normalization  t  this  peak  now  shows  some  of  the  sand 
reflections  appearing  as  clutter  on  the  image.  In  order  to  improve  the  image  before  a 
interpretation,  the  sand  reflection  was  removed.  The  interfering  sand  reflection  is 
illustrated  in  Figure  3.24  which  is  the  same  data  as  the  image  in  Figure  3.20  but  now 
normalized  to  the  peak  of  the  M19.  Thus  it  is  the  actual  surface  clutter  appearing  in  the 
M19  image.  If  the  threshold  is  now  raised  by  three  quantization  levels  so  the  lowest  level 


is  now  above  the  peak  surface  reflection,  an  improved  image  of  the  M19  now  results  as 
shown  in  Figure  3.25.  Now  the  detail  is  apparent  across  the  mine  surface  with  the  peak 
lobe  again  appearing  at  the  top  pressure  plate  (probably  due  to  the  air  gap),  with  slightly 
lower  lobes  at  the  top  sides  of  the  mine.  The  peak  again  occurs  at  the  actual  mine  depth. 

The  above  images  demonstrate  that  both  metal  and  plastic  buried  mines  can  be 
detected  and  shape  information  can  be  obtained  using  a  high  resolution  array  at  S-band  with 
bistatic  oblique  forward  scattring  near  the  Brewster  angle.  The  images  also  indicate  that  the 
signatures  of  the  two  mines  are  different  and  also  characteristic,  suggesting  the  ability  to 
identify  mine  type  by  the  image  signature. 

2.6.4  Conclusions 

The  results  of  this  task  are  very  important  since  they  demonstrate  the  feasibility  of 
detecting,  imaging  and  identifying  buried  metal  and  plastic  mines  using  a  focused  array  in  a 
bistatic  forward  scattering  geometry.  Further  improvement  in  the  images  could  be  obtained 
with  more  sophisticated  image  processing  techniques. 


Figure  3.20.  Focused  line  array  image  of  soil  reflections. 


transverse  distance 

Figure  3.21.  Focused  line  array  image  of  ruler  on  surface. 
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Figure  3.22.  Focused  line  array  image  of  Ml  5  at  7.5  cm  depth. 
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Figure  3.24.  Focused  line  array  image  of  soil  reflections  normalized  to  M19. 
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Figure  3.25.  Focused  line  array  image  of  M19  with  soil  reflections  removed. 
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2.7 


Propagation  Model 


2.7.1  Objective 

The  pupose  of  this  task  is  to  devleop  a  propagation  model  for  oblique  forward 
scattering  from  soil  and  to  compare  measured  data  with  the  model.  The  objective  is  to 
determine  the  effectiveness  of  the  Brewster  angle  effect  in  reducing  or  eliminating  soil 
surface  reflections. 

2.7.2  Propagation  Model 

A  theoretical  propagation  model  based  on  Stratton  [9]  has  been  developed  which 
predicts  soil  transmission  and  reflection  parameters  including  amplitude,  phase  and  index 
of  refraction  variations.  The  model  can  be  used  for  arbitrary  polarization,  angle  of 
incidence  and  soil  conductivity  and  dielectric  constant.  A  computer  model  has  been 
implemented  to  readily  compare  the  theory  with  experimental  measurements. 

Measurements  using  dry  and  moist  sand  have  been  performed  to  verity  the 
Brewster  angle  effect.  The  same  forward  scattering  geometry  is  used  as  for  the  previous 
experiments,  i.e.,  with  the  transmit  antenna  directed  forward  to  the  soil  suface  at  the  angle 
of  incidence,  with  the  receiver  symmetrically  located  at  the  angle  of  reflection  on  the  other 
side  of  the  incident  point.  The  angle  of  incidence  was  then  varied  for  the  measurements  by 
moving  the  transmitter  antenna  position  to  the  angle  of  reflection.  The  measurements  were 
performed  with  polarization  in  the  plane  of  incidence  (vertical)  since  a  vertically  polarized 
wave  at  the  Brewster  angle  would  be  completely  transmitted  into  the  soil. 

Figures  3.26  -  3.29  compare  the  results  of  the  theoretical  model  with  the  measured 
data.  They  show  the  reflection  coefficient  of  the  electric  field  as  a  function  of  the  angle  of 
incedence.  Figure  3.26  is  the  theoretical  result  for  L-band  with  sand  containing  5%  water. 
The  reflection  coefficient  becomes  very  small  between  60°  and  65°.  The  corresponding 
measured  data  is  shown  in  Figure  3.27  and  agrees  with  the  general  characteristics  of  the 
model.  Figure  3.28  shows  the  theoretical  S-band  result  for  dry  sand  having  zero  reflection 
coefficient  at  the  Brewster  angle  which  occurs  at  58°.  The  measured  data  of  Figure  3.29 
shows  good  agreement  and  has  a  minimum  reflection  coefficient  of  about  .03  at  55°. 

2.7.3  Conclusion 

This  task  has  demonstrated  conclusively  a  substantial  reduction  of  soil  surface 
reflections  for  a  foward  scattering  system  using  vertical  polarization  at  the  Brewster  angle. 
The  results  are  valid  for  both  dry  and  moist  soil,  the  moist  soil  having  a  minimum  reflection 
coefficient  slightly  above  zero  at  an  angle  somewhat  larger  than  the  Brewster  angle. 
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Figure  3.26.  Theoretical  soil  reflection  coefficient  at  L-band. 
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Figure  3.27.  Measured  soil  reflection  coefficient  at  L-band. 
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Figure  3.28.  Theoretical  soil  reflection  coefficient  at  S-band. 
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Figure  3.29.  Measured  soil  reflection  coefficient  at  S-band. 
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COMMON  TO  ALL  PROGRAMS 


(a)  Programing  Language 

Programs  have  been  written  in  FORTRAN  77  programing  language. 


(b)  Location 

The  programs  are  stored  in  the  account  named  ozdemir  on  the  Electrical  Engineering 
Department  branch  of  the  main  frame  computer  named  pender  located  at  the  School  of 
Engineering  and  Applied  Science  of  the  University  of  Pennsylvania.  They  are  also  available 
on  5"  and  3.5"  floppy  disks  provided. 


(c)  Execution 

(i)  To  execute  the  program  image./ ,  for  example,  on  the  main  frame  computer 
pender,  first,  to  compile  it,  type  "f77  image.f'.  If  no  errors  are  detected,  then  type  "a.out" 
to  run  it. 

(ii)  To  execute  a  program  on  an  IBM  PC,  first  move  the  program  file  and  all  the 
input  files  to  the  directory  thaUias  the  fortran  compiler  applications.  The  file  containing  the 
program  must  have  the  extention  ".for".  So,  to  execute  the  program  image.for,  for 
example,  type  "fori",  press  RETURN,  type  "image",  and  press  RETURN  four  times.  If 
no  errors  are  detected,  type  "pas2",  press  RETURN.  If  no  errors  are  detected,  type  "link", 
press  RETURN,  type  "image",  and  press  RETURN  four  times.  If  no  errors  are  detected, 
just  type  "image"  and  press  RETURN  to  run  the  program. 


(d)  Coordinate  Sign  Convention 


The  coordinate  sign  convention  is  given  in  Figure  1  and  used  by  the  programs 
described  here  for  identifying  the  air  and  the  subsurface  media.  . 


Figure  1:  Coordinate  Sign  Convension 


(e)  Flexibility 

The  programs  provide  quite  a  deal  of  flexibility  for  imaging  setup  in  the  sense  that 
transmitter,  scan  points  and  receiver  elements  can  be  in  either  medium  (air  or  subsurface 
medium.)  Only  serious  constraints  upon  the  imaging  setup  are  that  the  transmitter  is  a 
single  point  element,  and  that  the  receive  array  is  linear,  periodic  and  parallel  to  one  of  x,  y 
or  z  axes. 

The  way  they  are  presented  here,  the  programs  are  consistent  with  each  other  in 
terms  of  formating  for  the  flow  of  data  between  the  programs.  The  samples  of  input  and 
output  files  are  available  in  the  respective  directories.  The  sample  data  in  the  input  files  are 
in  accordance  with  the  current  imaging  setup  and  the  system  parameters. 


(f)  Units 


MKS  unit  system  is  used.  Any  exceptions  are  cited  wherever  appear. 
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PROGRAM  #1:  interface ./ 


Location :  ozdemir  I  research! interface 

1.1  CODE 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


THIS  SUBROUTINE  CALCULATES  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE 
PARTS  OF  THE  RAY  PATH  BETWEEN  TWO  SEPERATE  POINTS  (PI  AND  P2)  IN  A  MEDIA 
COMPOSED  OF  AIR  AND  A  SUBSURFACE  MEDIUM,  EACH  OF  WHICH  IS  A  HOMOGENOUS, 
ISOTROPIC,  NONMAGNETI ZABLE  AND  CONDUCTING  (EXCEPT  FOR  AIR,  WHICH  IS 
NONCONDUCTING)  MEDIUM,  SEPERATED  BY  A  PLANAR  INTERFACE.  THE  POINTS  CAN  BE 
IN  DIFFERENT  MEDIUMS  AS  WELL  AS  IN  THE  SAME  MEDIUM.  TEE  INTERFACE  IS  THE 
Z=0  OR  (X-Y)  PLANE.  POINTS  WITH  Z<=0  ARE  IN  THE  AIR  MEDIUM  WHILE  POINTS 
WITH  Z>0  ARE  IN  THE  SUBSURFACE  MEDIUM. 


THE  PARAMETERS  AND  THE  VARIABLES  USED  IN  THE  PROGRAM: 


eps 

(INPUT 

PARAMETER)  : 

sig 

( 

It 

) 

f 

( 

II 

) 

xl 

( 

It 

) 

yi 

( 

It 

) 

zl 

( 

II 

) 

x2 

( 

11 

) 

y2 

( 

11 

) 

z2 

( 

11 

) 

nit 

( 

II 

) 

eta 

(OUTPUT  PARAMETER) : 

dl 

( 

11 

)  : 

d2 

( 

11 

)  : 

RELATIVE  DIELECTRIC  CONSTANT  OF  THE  SUBSURFACE 
MEDIUM  (>=1) 

CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter] 
FREQUENCY  OF  RADIATION  [GHz] 

X-COORDINATE  OF  THE  FIRST  POINT  (PI)  [cm] 

Y  " 

Z  « 

X-COORDINATE  OF  THE  SECOND  POINT  (P2)  [cm] 

Y  " 

Z  " 

NUMBER  OF  ITERATIONS  IN  CALCULATING  THE  INCIDENCE 
ANGLE 

REFRACTIVE  INDEX  OF  THE  SUBSURFACE  MEDIUM 
EVALUATED  AT  THE  CALCULATED  INCIDENCE  ANGLE 
■ME  LENGTH  OF  THE  AIR  PART  OF  THE  RAY  PATH  [cm] 
THE  LENGTH  OF  THE  SUBSURFACE  PART  OF  THE  RAY  PATH 
[cm] 


theta  (INTERMEDIATE  VARIABLE)  :  INCIDENCE  ANGLE  [rad] 
phi  (  "  )  :  REFRACTION  ANGLE  [rad] 

d  (INTERMEDIATE  PARAMETER) :  HORIZONTAL  DISTANCE  EETWEEN  THE  POINTS  PI 

AND  P2  [cm] 

h  (  "  ):  HIGHT  OF  THE  POINT  (Fl  OR  P2)  IN  THE  AIR 

FROM  THE  INTERFACE  [cm] 

z  (  "  )  :  DEPTH  OF  THE  POINT  IN  THE  SUBSURFACE 

MEDIUM  FROM  THE  SURFACE  (INTERFACE)  [cm] 


subroutine  interface (eps, sig, f , xl, yl,  zl ,  x2,  y2 ,  z2,  nit ,  eta, dl, d2) 
real  eps, sig, f , xl, yl, zl, x2, y2, z2, eta, dl,d2 
integer  nit 

REDUCTION  OF  THE  3-D  PROBLEM  TO  A  PLANAR  GEOMETRY  PROBLEM 
d=sqrt ( (xl-x2)**2.+(yl-y2)**2.) 

IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  AIR  MEDIUM  OR  eps=l 

if((zl.le.O.. and. z2 . le . 0 . ) .or .eps .eq. 1 . )  then 
dl=sqrt (d**2 .  +  (zl-z2) **2  . ) 

d2=0 . 


eta  =  l . 

goto  3  5 

*  '  end  if 

c 

c  IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  SUBSURFACE  MEDIUM.  .  . 

c 

*  if ( zl . ge . 0 . . and . z2 . ge . 0 . )  then 
dl=0 . 

d2=sqrt (d**2 . + (zl-z2) **2 . ) 

eta=l . /sqrt (2 . ) *sqrt (eps+sqrt (eps**2 .+(18.*sig/f)*2.)) 
goto  3 
end  if 
c 

C  IF  THE  POINTS  ARE  IN  DIFFERENT  MEDIA 

C 

if (zl . le . 0 . .and. z2 . gt . 0 . )  then 
h=-zl 
z  =  z2 
goto  1 
else 
h=-z2 
z=zl 
end  if 
c 

C  ITERATION  FOR  THE  INCIDENCE  ANGLE  BEGINS  WITH  AN  INITIAL  GUESS 

c 

1  theta=atan (d/ (h+z/eps ) ) 
do  2  it=l,nit 

eta=l . /sqrt (2 . ) *sqrt (eps+sin (theta) **2 . +sqrt ( (eps- 
%  sin (theta) **2.)**2.  +  (18.  *sig/f ) **2 . ) ) 
phi=asin (sin (theta) /eta) 

2  theta=atan ( (d-z*tan (phi) ) /h) 
c 

C  "nit"  th  ITERATION  VALUES  FOR  THE  INCIDENCE  AND  THE  REFRACTION  ANGLES 

c  "theta"  AND  "phi"  AND  THE  REFRACTIVE  INDEX  "eta"  ARE  TAKEN  TO  BE  THEIR 

C  SOLUTIONS .THEN  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  PARTS  OF  THE  RAY 

c  PATH  FOLLOW 

c 

dl=h/cos (theta) 
d2=z/cos (phi) 
return 
end 


3 
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1.2  GENERAL  INSTRUCTIONS  AND  COMMENTS 

This  is  not  a  directly  executable  routine  but  a  subroutine  to  be  appended  to  the  end 
of  the  main  programs  and  called  appropriately.  It  performs  the  most  basic  calculations 
required  by  our  wave  propagation  model  [1]. 

1.3  DESCRIPTION  OF  THE  PROBLEM  AND  THE  METHODS  OF 

SOLUTION  USED  BY  THE  SUBROUTINE 

The  most  basic  calculation  required  by  our  wave  propagation  model  is  determining 
the  lengths  of  the  air  and  the  subsurface  parts  of  the  total  ray  path  between  two  given  points 
(Pi  and  P2)  located  anywhere  in  a  media  composed  of  air  and  a  subsurface  medium, 
separated  by  a  planar  interface  (each  medium  is  homogeneous,  isotropic,  nonmagnetizable 
and  conducting  except  for  air,  which  is  nonconducting,)  i.e.,  calculating  di  and  d2  shown 
in  Figure  2. 


Figure  2:  3-D  Ray  Propagation  Geometry 

Although,  at  first,  the  problem  looks  like  a  3-D  problem,  the  propagation  of  the  ray 
from  Pi  to  P2  takes  place  in  the  incidence  plane  defined  by  the  air  and  the  subsurface  parts 


of  the  total  ray  path,  and  therefore,  it  is  actually  a  2-D  problem.  Then  consider  the  planar 
geometry  in  Figure  3,  where,  with  reference  to  Figure  2, 


h  =  -Zj 

z  =  z2 

d =V  (x2~  xi)2+  (y2-yi)2 


From  the  Snell's  Law,  we  have 


where 


Sin  9  = 


Sin  8 
*1(0) 


Sin20  +  yj (er-  Sin2e)2+ 


(1) 


(2) 


and  is  called  the  refractive  index  of  the  subsurface  medium.  £j  and  c  are  the  relative 
dielectric  constant  and  the  conductivity  of  the  subsurface  medium,  respectively,  and  f  is 
the  frequency  of  radiation  in  units  of  GHz. 

From  the  geometry  in  Figure  3, 
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d  =  h  tanO  +  z  tantp 


Substituting  (1)  in  (3), 


[Sin9 

d  =  h  tan 9  +  z  tan  Sin*1 


where  the  only  unknown  is  the  incidence  angle  0.  The  difficulty  of  obtaining  an  analytical 
solution  for  0  from  the  last  equation  is  surely  appreciated.  Then,  we  have  searched  for  an 
iterative  solution  and  come  up  with  the  following  solution,  which  has  been  observed  to 
converge  quite  fast:  We  start  with  a  good  initial  guess,  which  we  have  determined,  by 
trying,  to  be 


0(O)  =  tan-^h  +dz/er 


Then  using  (1)  and  (3),  we  have  the  following  iterations  for  cp  and  0: 


<P(n)  =  Sin-i 


SinQfr-1) 

ti[0(n-i>] 


,  .  ,["d-ztan(p(n> 

6(n)=  tan-i[ - -J  ,  n  =  1,2,3,... 


where  ri(0)  is  as  given  in  (2). 


For  z<50  cm,  0  converges  reasonably  close  to  the  solution  in  less  than  5  iterations. 
Having  calculated  9  and  6>  di  and  6.2  can  be,  from  Figure  3,  simply  expressed  as 


Cos0 


d2  =  -^- 
2  Coscp 


Also  the  refractive  index  tj(0)  is  readily  computed  from  (2). 
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PROGRAM  #2:  proptime.f 
Location :  ozdemir! research! proptime 

2.1  CODE 


c  This  program  calculates  the  total  propagation  time  of  the  wave  from  the 

c  transmitter  to  each  specified  scan  point  and  then  from  each  scan  point  to 

c  each  specified  receiver  element  in  a  media  composed  of  air  and  a 

c  conductive  subsurface  medium  seperated  by  a  planar  interface.  There  is 

c  only  one  transmitter  element.  Receiver  array  is  periodic  and  oriented 

c  parallel  to  one  of  x,y  or  z  axes, 

c 

c  PARAMETERS  AND  VARIABLES  USED  IN  THE  PROGRAM 

c 

c  INPUT  PARAMETERS: 

c 

c  eps : 

c  sig : 

c  f : 

c  nit: 

c 

c  nel : 

c  spel: 

c  xrc : 

c  yrc : 

c  zrc : 

c  dir : 

c  xtr: 

c  ytr : 

c  ztr : 

c  zlsc: 

c  z2sc: 

c  dzsc: 

c  ylsc: 

c  y2sc: 

c  dysc: 

c  xlsc: 

c  x2sc: 

c  dxsc : 

c 
c 

c  fname: 

c 


c 


c 

VARIABLES : 

c 

c 

nzsc: 

Number  of  scan  points  in  z 

direction 

c 

nysc: 

II 

y 

It 

c 

nxsc: 

tl 

z 

II 

c 

zsc : 

z-coordinate  of 

the  current 

scan  point  [cm] 

c 

ysc : 

y 

tl 

c 

xsc : 

X 

II 

c 

xel : 

x-coordinate  of 

the  current 

receiver  element  position  [cm] 

c 

yel: 

y 

II 

c 

zel : 

z 

II 

c 

del : 

Distance  of  the 

current  receiver  element  position  to  the  center 

c 

of  the  receiver 

a  r  r  a  v  f  cm  1 

Relative  dielectric  constant  of  the  subsurface  medium  (>=1) 
Conductivity  of  the  subsurface  medium  [mho/meter] 

Frequency  of  radiation  [GHz] 

Number  of  iterations  for  subroutine  "interface"  to  calculate  the 
incidence  angle 

Number  of  elements  of  the  receiver  array 

Constant  receiver  element  spacing  [cm] 

x-coordinate  of  the  center  of  the  receiver  array  [cm] 


direction  (orientation)  of  the  receiver  array  (x,y  or  z) 
x-coordinate  of  the  transmitter  [cm] 


z-coordinate  of  the  first  scan  point  [cm] 

"  ibst  " 

Spacing  between  the  scan  points  in  z  direction  [cm]  (non  zero) 
y-coordinate  of  the  first  scan  point  [cm] 

"  last 

Spacing  between  the  scan  points  in  y  direction  [cm]  (non  zero) 
x-coordinate  of  the  first  scan  point  [cm] 

"  last  " 

Spacing  between  the  scan  points  in  x  direction  [cm]  (non  zero) 

(If  no  scanning  is  desired  in,  for  example,  x  direction,  specify  in 
the  input  file  xlsc=x2sc  and  dxsc=l) 

Name  of  the  output  file  for  storing  the  computed  propagation  times 
(maximum  25  characters) 


c 

c 

<  c 

c 

c 

c 

*  c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Refractive  index  of  the  subsurface  medium  evaluated  at  the 
incidence  angle  satisfying  the  ray  geometry  in  transmit  or  receivejQ 
incidence  plane 

Length  of  the  air  part  of  the  ray  path  in  transmit  or  receive 
incidence  plane 

Length  of  the  subsurface  part  of  the  ray  path  in  transmit  or 
receive  incidence  plane 

Propagation  time  in  the  transmit  incidence  plane  [nsec] 

"  receive  " 

Total  propagation  time  (tl  +  t2)  [nsec] 


dimension  t ( 5 0 ) 
character *2 5  fname 
character  dir 

OPEN  THE  INPUT  FILE  "proptime .  in"  AND  READ  THE  INPUT  PARAMETERS 
DESCRIBING  THE  IMAGING  SET-UP  AND  THE  MEDIA 

open  (11, file=' proptime . in' ) 
read(ll, '  (flO .4) ' )  eps,sig,f 
read(ll, '  (i3) ' )  nit,nel 
read(ll,'  (fl0.4)')  spel, xrc, yrc, zrc 
read (11, ' (a) ' )  dir 

read(ll,'  (fl0.4)')  xtr,  ytr, ztr, zlsc, z2sc, dzsc, ylsc, y2sc,  dysc, 

%  xlsc, x2sc, dxsc 
read (11, ' (a25) ' )  fname 

DETERMINING  THE  NUMBER  OF  SCAN  POINTS  IN  z,y  AND  x  DIRECTIONS 

nzsc=nint ( (z2sc-zlsc) /dzsc) +1 
nysc=nint ( (y2sc-ylsc) /dysc) +1 
nxsc=nint ( (x2sc-xlsc) /dxsc) +1 

OPEN  THE  OUTPUT  FILE  FOR  STORING  THE  PROPAGATION  TIMES 
open (21, f ile=fname, status=' new' ) 

SPECIFYING  THE  LOCATIONS  OF  THE  SCAN  POINTS  IN  x,  y  and  z  DIRECTIONS 

do  1  ixsc=l,nxsc 

xsc=xlsc+dxsc*float (ixsc-1) 
do  1  izsc=l,nzsc 

zsc=zlsc+dzsc*f loat (izsc-1) 
do  1  iysc=l,nysc 

ysc=ylsc+dzsc*f loat (iysc-1) 

CALCULATING  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  PARTS  OF  THE  RAY 
PATH  FROM  THE  TRANSMITTER  TO  THE  SCAN  POINT 

call  interface (eps, sig, f ,  xtr, ytr , ztr, xsc, ysc, zsc, nit, eta,  dl ,  d2) 
CALCULATING  THE  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE  SCAN  POINT 
tl=(dl+eta*d2)/30. 

SPECIFYING  RECEIVER  ELEMENT  POSITIONS 

xel=xrc 

yel=yrc 

zel=zrc 

do  2  irc=l,nel 

del=spel* (float (ire) -float (nel+1) /2 . ) 
if (dir.eq.'x' )  xel=xrc+del 

'f(di>'.pa,'vM  v<=l 


eta : 

dl: 

d2 : 

tl: 
t2 : 
t : 


if (dir . eq . ' z ' )  zel=zrc+del 
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* 


c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

2 

c 

c 

c 

c 

1 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CALCULATING  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  FARTS  OF  THE  RAY 
PATH  FROM  THE  SCAN  POINT  TO  THE  RECEIVER  ELEMENT 

call  interface (eps, sig,  f ,  xsc, ysc, zsc, xel, yel,  zel,  nit,  eta, dl, d2) 

CALCULATING  THE  PROPAGATION  TIME  FROM  THE  SCAN  POINT  TO  THE  RECEIVER 

t2=(dl+eta*d2)/30. 

CALCULATING  THE  TOTAL  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE 
SCAN  POINT  AND  FROM  THERE  TO  THE  RECEIVER 

t  (irc)=tl+t2 

WRITE  INTO  THE  OUTPUT  FILE  THE  TOTAL  PROPAGATION  TIMES  FOR  THE  CURRENT 
SCAN  POINT 

write  (21,'  (8f 9 . 5) ' )  (t(irc) , irc=l, nel) 

stop 

end 

THIS  SUBROUTINE  CALCULATES  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE 
PARTS  OF  THE  RAY  PATH  BETWEEN  TWO  SEPERATE  POINTS  (PI  AND  P2)  IN  A  MEDIA 
COMPOSED  OF  AIR  AND  A  SUBSURFACE  MEDIUM,  EACH  OF  WHICE  IS  A  HOMOGENOUS, 
ISOTROPIC,  NONMAGNETI ZABLE  AND  CONDUCTING  (EXCEPT  FOR  AIR,  WHICH  IS 
NONCONDUCTING)  MEDIUM,  SEPERATED  BY  A  PLANAR  INTERFACE.  THE  POINTS  CAN  BE 
IN  DIFFERENT  MEDIUMS  AS  WELL  AS  IN  THE  SAME  MEDIUM.  TEE  INTERFACE  IS  THE 
Z=0  OR  (X-Y)  PLANE.  POINTS  WITH  Z<=0  ARE  IN  THE  AIR  MEDIUM  WHILE  POINTS 
WITH  Z>0  ARE  IN  THE  SUBSURFACE  MEDIUM. 


THE  PARAMETERS  AND  THE  VARIABLES  USED  IN  THE  CODE: 

eps  (INPUT  PARAMETER)  :  RELATIVE  DIELECTRIC  CONSTANT  CF  THE  SUBSURFACE 

MEDIUM  (>=1) 

CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter] 
FREQUENCY  OF  RADIATION  [GHz] 

X-COORDINATE  OF  THE  FIRST  POINT  (PI)  [cm] 

Y—  " 

Z  " 

X-COORDINATE  OF  THE  SECOND  POINT  (P2)  [cm] 

Y  " 

Z  " 

NUMBER  OF  ITERATIONS  IN  CALCULATING  THE  INCIDENCE 
ANGLE 

eta  (OUTPUT  PARAMETER)  :  REFRACTIVE  INDEX  OF  THE  SUBSUFE’ACE  MEDIUM 

EVALUATED  AT  THE  CALCULATED  INCIDENCE  ANGLE 
dl  (  "  )  :  THE  LENGTH  OF  THE  AIR  PART  OF  THE  RAY  PATH  [cm] 

d2  (  "  )  :  THE  LENGTH  OF  THE  SUBSURFACE  FART  OF  THE  RAY  PATH 

[cm] 

theta  (INTERMEDIATE  VARIABLE)  :  INCIDENCE  ANGLE  [rad] 
phi  (  "  )  :  REFRACTION  ANGLE  [rad I 

d  (INTERMEDIATE  PARAMETER) :  HORIZONTAL  DISTANCE  BETWEEN  THE  POINTS  PI 

AND  P2  [cm] 

h  (  "  ) :  HIGHT  OF  THE  POINT  (PI  OR  P2)  IN  THE  AIR 

FROM  THE  INTERFACE  [cm] 

z  (  "  )  :  DEPTH  OF  THE  POINT  IN  THE  SUBSURFACE 

MEDIUM  FROM  THE  SURFACE  (INTERFACE)  [cm] 


sig  v 

( 


f 

xl 

yi 

zl 
x2 
y2 

z2  ( 
nit  ( 


subroutine  interface (eps,  sig, f , xl,yl, zl,x2,  y2,  z2,  nit,  eta,  dl,  d2) 
real  eps, sig, f , xl, yl,  zl, x2, y2, z2, eta, dl,d2 
integer  nit 


c 


12 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


1 


2 


c 

c 

c 

c 

c 

c 


REDUCTION  OF  THE  3-D  PROBLEM  TO  A  PLANAR  GEOMETRY  PROBLEM 
d=sqrt ( (xl-x2)**2.+(yl-y2)**2.) 

IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  AIR  MEDIUM  OR  eps=l..., 

if ( (zl . le . 0 . . and. z2 . le . 0 . ) . or .eps . eq. 1 . )  then 
dl=sqrt (d**2 .+ (zl-z2) **2. ) 
d2=0 . 
eta=l . 
goto  3 
end  if 

IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  SUBSURFACE  MEDIUM... 

if (zl .ge . 0 . .and. z2 .ge . 0 . )  then 
dl=0 . 

d2=sqrt (d**2 .+ (zl-z2) **2 . ) 

eta=l . /sqrt (2 . ) *sqrt (eps+sqrt (eps**2 .+ (18 . *sig/f ) *2 . ) ) 
goto  3 
end  if 

IF  THE  POINTS  ARE  IN  DIFFERENT  MEDIA 

if (zl . le . 0 . .and. z2 . gt . 0 . )  then 
h=-zl 
z=z2 
goto  1 
else 
h=-z2 
z=zl 
end  if 

ITERATION  FOR  THE  INCIDENCE  ANGLE  BEGINS  WITH  AN  INITIAL  GUESS 

theta=atan (d/ (h+z/eps) ) 
do  2  it*=l,nit 

eta=l . /sqrt (2 . ) *sqrt (eps+sin (theta) **2 . +sqrt ( (eps- 
%  sin (theta) **2. ) **2 .  + (18 . *sig/f ) **2 . ) ) 
phi=asin (sin (theta) /eta) 
theta=atan ( (d-z*tan (phi)  )Vh) 

"nit"  th  ITERATION  VALUES  FOR  THE  INCIDENCE  AND  THE  REFRACTION  ANGLES 
"theta"  AND  "phi"  AND  THE  REFRACTIVE  INDEX  "eta"  ARE  TAKEN  TO  BE  THEIR 
SOLUTIONS  .THEN  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  PARTS  OF  THE  RAY 
PATH  FOLLOW 

dl=h/cos (theta) 
d2=z/cos (phi) 
return 
end 
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2.2  GENERAL  INSTRUCTIONS  AND  COMMENTS 


The  program  makes  use  of  the  previously  described  subroutine  interface ./  (in 
calculating  the  lengths  of  the  refracted  ray  paths),  which  has  been  appended  to  the  end  of 
the  program.  It  reads  the  input  data  from  an  input  file  called  proptime. in. 

To  run  this  program  on  IBM  PC,  just  change  the  extention  of  the  program  file  name 
from  "f  to  "for",  and  follow  the  instructions  given  in  the  beginning  of  this  report  on 
running  programs  on  IBM  PCs. 

The  runing  time  of  the  program  is  roughly  proportional  to  the  number  of  scan 
points  times  the  number  of  receiver  elements.  Let  the  number  of  scan  points  in  longitudinal, 
vertical  and  horizontal  directions  be  Nx,  Nz,  Ny,  and  the  number  of  receiver  elements  be 
Nr,  respectively.  When  Nx=l,  Nz=51,  Ny=101,  Nr=22,  the  program  took  126  minutes  to 
run.  Then  the  formula  for  runing  time  (on  IBM  PC)  is 

Runing  Time  =  126  x(NxxNzxNyxNr)/(  1  x51  x  101  x  22  )  Minutes. 


2.2.1  INPUT  FILE:  proptime.in 


Format: 


Input 

mf 

LingJ  Parameter  Egmiai 

Description  of  Parameter 

1 

eps 

fl0.4 

Relative  Dielectric  Constant  of  the  Subsurface  Medium  (>1) 

2 

sig 

t» 

Conductivity  of  the  Subsurface  Medium  [mho/meter]  (>0) 

3 

f 

ti 

Frequency  of  Radiation  [GHz] 

4 

nit 

i3 

Number  of  Iterations  for  Subroutine  "interface"  to  Calculate  the 

Incidence  Angle 

5 

nel 

it 

Number  of  Elements  of  the  Receiver  Array 

6 

spel 

fl0.4 

Constant  Receiver  Element  Spacing  [cm] 

7 

xrc 

if 

x-Coordinate  of  the  Center  of  the  Receiver  Array  [cm] 

8 

yrc 

tt 

y 

9 

zrc 

it 

z 

10 

dir 

a 

Direction  of  the  Receiver  Array:  Put  Y  If  Parallel  to  x 

Put  y  "  y 


Put  'z' 


z 


11 

xtr 

fl0.4 

x-Coordinate  of  the  Transmitter  [cm] 

12 

ytr 

11 

y 

13 

ztr 

it 

z 

14 

zlsc 

11 

z-Coordinate  of  the  First  Scan  Point  [cm] 

15 

z2sc 

M 

Last 

16 

dzsc 

11 

Spacing  Between  the  Scan  Points  in  z  Direction  [cm]  (* 0 ) 

17 

ylsc 

H 

y-Coordinate  of  the  First  Scan  Point  [cm] 

18 

y2sc 

11 

"  Last 

19 

dysc 

it 

Spacing  Between  the  Scan  Points  in  y  Direction  [cm]  (5O) 

20 

xlsc 

it 

x-Coordinate  of  the  First  Scan  Point  [cm] 

21 

x2sc 

11 

"  Last 

22 

dxsc 

11 

Spacing  Between  the  Scan  Points  in  x  Direction  [cm]  0*0) 

(If  no  scaning  (imaging)  is  desired  in,  for  example,  x 
direction,  specify  in  the  input  file  xlsc=x2sc  and  dxsc=lcm) 

23 

fname 

a25 

Name  of  the  Output  File  for  Storing  the  Computed 

Propagation  Times  (maximum  25  characters) 


2.3  DESCRIPTION  OF  THE  PROBLEM  AND  THE  METHODS  OF 
SOLUTION  USED  BY  THE  PROGRAM 

This  program  calculates,.the  total  propagation  time  of  the  wave  from  the  transmitter 
location  to  each  specified  scan  point  and  then  from  the  scan  point  to  each  specified  receiver 
position  in  a  media  composed  of  air  and  a  conductive  subsurface  medium  separated  by  a 
planar  interface.  Transmitter  location  and  the  scan  points  as  well  as  the  receiver  positions 
can  be  specified  in  any  medium  (air  or  subsurface  medium.)  The  restrictions  are:  The 
transmitter  has  one  element,  the  scan  points  can  span  a  3-D  space,  and  the  receiver 
positions  should  be  lined  up  parallel  to  one  of  x,  y  or  z-axes.  Consider  the  geometry  in 
Figure  4. 

As  it  is  seen  in  Figure  4,  the  total  propagation  time  consists  of  two  parts.  One  in  the 
transmit  incidence  plane  and  the  other  in  the  receive  incidence  plane.  The  geometry  of  the 
problem  in  each  plane  is  basically  the  same  as  shown  in  Figure  2,  where  the  problem  was 
to  compute  the  lengths  of  the  air  and  the  subsurface  parts  of  the  ray  path,  i.e.,  di  and  62, 
which  is  exactly  what  the  subroutine  interface ./  does.  Then  we  can  use  this  subroutine  to 


compute  di  and  d2  in  both  transmit  and  receive  incidence  planes,  and  with  the  knowledge 
of  the  propagation  speeds  in  each  medium,  which  are  determined  by  the  relative  dielectric 
constant  (£r)  and  the  conductivity  (a)  of  the  respective  medium  through  the  refractive  index, 
we  can  compute  the  total  propagation  time. 


Figure  4:  3-D  Ray  Propagation  Geometry 


Let  the  lengths  of  the  air  and  the  subsurface  parts  of  the  ray  path  in  the  Transmit  and 
Receive  Incidence  Planes  be  duand  d2t,  and  djr  and  d2r,  and  the  subsurface  propagation 
speeds  be  c/ri(0t)  and  c/r|(8r),  respectively,  c  is  the  speed  of  light  in  free-space.  Ti(0t)  and 
T|(0r)  are  the  refractive  indices  of  the  subsurface  medium  evaluated  at  the  incidence  angles 
satisfying  the  ray  geometry  in  Transmit  and  Receive  Incidence  Planes,  respectively.  Then 
the  total  propagation  time  from  the  transmitter  to  the  receiver  is 


16 


t 


d2t  ,  d2r  ,  dir 


Tl(0t)  Tl(0r) 


(6) 


This  is  repeated  for  each  scan  point  and  receiver  position  specified 

As  long  as  the  parameters  and  the  geometry  of  the  imaging  setup  remain  constant, 
the  propagation  times  will  not  change.  Therefore,  they  need  to  be  computed  only  once. 
After  they  are  stored,  they  can  be  reused  for  imaging  of  different  targets  (by  the  program  to 
be  explained  next)  for  which  the  I  and  Q  data  are  collected  with  the  same  imaging  setup. 
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PROGRAM  #3:  image.f 
Location:  ozdemir!  research! imaging 


3.1  CODE 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  PROGRAM  FORMS  THE  IMAGE  OF  A  SPECIFIED  TARGET  REGION  USING  THE 
P RE -CALCULATED  PROPAGATION  TIMES  AND  THE  TIME  SAMPLES  OF  THE  INPHASE  AND 
THE  QUADRATURE  COMPONENTS  OF  THE  COMPLEX  WAVEFORMS  MEASURED  AT  THE 
OUTPUT  OF  THE  RECEIVER  ARRAY  ELEMENTS.  THE  IMAGE  FORMING  IS  PERFORMED  IN 
ACCORDANCE  WITH  THE  IMAGING  MODEL  PRESENTED  IN  "VFRC  PROGRESS  REPORT  NO. 
59"  (REFERENCE  1) . 


INPUT  PARAMETERS: 


nel :  NUMBER  OF  RECEIVER  ARRAY  ELEMENTS 

nsamp :  NUMBER  OF  TIME  SAMPLES  OF  THE  I  AND  Q  WAVEFORMS 
nz:  NUMBER  OF  SCAN  POINTS  IN  z-DIRECTION 


ny : 
nx : 


It 

It 


y 

X 


tt 

It 


iref:  THE  SAMPLE  NUMBER  CORRESPONDING  TO  THE  CENTER  OF  THE  PULSE 

SCATTERED  FROM  THE  REFERENCE  SCATTERING  POINT  AND  MEASURED  AT  THE 
REFERENCE  MEASUREMENT  POINT, 
xtr:  x-COORDINATE  OF  THE  TRANSMITTER  [cm] 
ytr :  y  " 

ztr:  z  " 


xrs : 
yrs : 
zrs : 
xrm: 
yrm: 
zrm: 


x-COORDINATE  OF  THE  REFERENCE  SCATTERING  POINT  [cm] 

y  " 

z  " 

X-COORDINATE  OF  THE  REFERENCE  MEASUREMENT  POINT  [cm] 


eps:  RELATIVE  DIELECTRIC  CONSTANT  OF  THE  SUBSURFACE  MEDIUM  (>=1) 
sig:  CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter]  (>=0) 
f:  FREQUENCY  OF  RADIATION  [GHz] 

tsamp:  SAMPLING  TIME  INTERVAL  FOR  I  AND  Q  DATA  [nsec] 
iqname:  NAME  OF  THE  INPUT  FILE  CONTAINING  THE  MEASURED  I  AND  Q  DATA 
tmname :  NAME  OF  THE  INPUT  FILE  CONTAINING  THE  COMPUTED  PROPAGATION  TIMES 
imname :  NAME  OF  THE  OUTPUT  FILE  FOR  STORING  THE  PRODUCED  IMAGE  DATA 


VARIABLES : 

in:  INPHASE  COMPONENT  OF  THE  RECEIVED  FIELD 
qd:  QUADRATURE  COMPONENT  OF  THE  RECEIVED  FIELD 

tl:  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE  REFERENCE  SCATTERER 
POINT  [nsec] 

t2 :  PROPAGATION  TIME  FROM  THE  REFERENCE  SCATTERER  POINT  TO  THE  REFERENCE 
MEASUREMENT  POINT  [nsec] 

tref:  REFERENCE  TOTAL  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE 
REFERENCE  MEASUREMENT  POINT  (tl  +  t2)  [nsec] 
w:  ANGULAR  FREQUENCY  OF  RADIATION  [rad/nsec] 

t:  TOTAL  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE  SPECIFIED  SCAN 

POINT  AND  THEN  TO  THE  SPECIFIED  RECEIVE  ARRAY  ELEMENT  POSITION  [nsec] 
isamp:  SAMPLE  NUMBER  CORRESPONDING  TO  THE  CENTER  OF  THE  PULSE  SCATTERED 
FROM  THE  SPECIFIED  SCAN  POINT  AND  MEASURED  AT  THE  SPECIFIED 
RFPFTVE  ARRAY  FJ.FMFNT 


c 


*  C 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


s (nx, nz , ny) =150 . 

OPEN  THE  OUTPUT  FILE  AND  STORE  THE  IMAGE  DATA 
open (51 , f ile=imname) 

write <51, '  (10li4) ' )  (((nint(s(ix,iz,iy)),  iy=l,  ny) , iz=l, nz) , 
%  ix=l,nx) 
stop 
end 


THIS  SUBROUTINE  CALCULATES  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE 
PARTS  OF  THE  RAY  PATH  BETWEEN  TWO  SEPERATE  POINTS  (PI  AND  P2)  IN  A  MEDIA 
COMPOSED  OF  AIR  AND  A  SUBSURFACE  MEDIUM,  EACH  OF  WHICH  IS  A  HOMOGENOUS, 
ISOTROPIC,  NONMAGNETI ZABLE  AND  CONDUCTING  (EXCEPT  FOR  AIR,  WHICH  IS 
NONCONDUCTING)  MEDIUM,  SEPERATED  BY  A  PLANAR  INTERFACE.  THE  POINTS  CAN  BE 
IN  DIFFERENT  MEDIUMS  AS  WELL  AS  IN  THE  SAME  MEDIUM.  THE  INTERFACE  IS  THE 
Z=0  OR  (X-Y)  PLANE.  POINTS  WITH  Z<=0  ARE  IN  THE  AIR  MEDIUM  WHILE  POINTS 
WITH  Z>0  ARE  IN  THE  SUBSURFACE  MEDIUM. 


THE  PARAMETERS  AND  THE  VARIABLES  USED  IN  THE  CODE: 


eps  (INPUT  PARAMETER) 

sig  (  "  ) 

f  (  ) 

xl  (  *'  ) 

yl  (  "  ) 

zl  (  •'  ) 

x2  (  •'  ) 

y2  (  "  ) 

Z2  (  "  ) 

nit  (  "  ) 

eta  (OUTPUT  PARAMETER) 


dl 

d2 


theta 

phi 

d 


) 


(INTERMEDIATE  VARIABLE) 

(  "  ) 
(INTERMEDIATE  PARAMETER) 


( 


RELATIVE  DIELECTRIC  CONSTANT  OF  THE  SUBSURFACE 
MEDIUM  (>=1) 

CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter] 
FREQUENCY  OF  RADIATION  [GHz] 

X-COORDINATE  OF  THE  FIRST  POINT  (PI)  [cm] 

Y  " 

Z  " 

X-COORDINATE  OF  THE  SECOND  POINT  (P2)  [cm] 

Y  " 

Z  " 

NUMBER  OF  ITERATIONS  IN  CALCULATING  THE  INCIDENCE 
ANGLE 

REFRACTIVE  INDEX  OF  THE  SUBSURFACE  MEDIUM 
EVALUATED  AT  THE  CALCULATED  INCIDENCE  ANGLE 
THE  LENGTH  OF  THE  AIR  PART  OF  THE  RAY  PATH  [cm] 
THE  LENGTH  OF  THE  SUBSURFACE  PART  OF  THE  RAY  PATH 
[etn] 

INCIDENCE  ANGLE  [rad] 

REFRACTION  ANGLE  [rad] 

HORIZONTAL  DISTANCE  BETWEEN  THE  POINTS  PI 
AND  P2  [cm] 

HIGHT  OF  THE  POINT  (PI  OR  P2)  IN  THE  AIR 
FROM  THE  INTERFACE  [cm] 

DEPTH  OF  THE  POINT  IN  THE  SUBSURFACE 
MEDIUM  FROM  THE  SURFACE  (INTERFACE)  [cm] 


subroutine  interface (eps, sig, f , xl, yl, zl, x2, y2, z2, nit,  eta,  dl,  d2) 
real  eps, sig,f,xl,yl,zl,x2,y2,z2,eta,dl,d2 
integer  nit 

REDUCTION  OF  THE  3-D  PROBLEM  TO  A  PLANAR  GEOMETRY  PROBLEM 
d=sqrt ( (xl-x2) **2 .  +  (yl-y2) **2  . ) 

IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  AIR  MEDIUM  OR  eps=l 

if ( (zl . le . 0 . .and. z2.1e.O.) .or. eps . eq. 1 . )  then 
dl=sqrt (d**2.+(zl-z2) **2. ) 
d2=0 . 
eta=l . 

CTOt  o  3 


o  o  o  o 


c 

c 

c 


end  if 


20 


i 


c 


IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  SUBSURFACE  MEDIUM... 

if  (zl .ge . 0 . .and . z2 . ge . 0 . )  then 
dl=0 . 

d2=sqrt (d**2.+(zl-z2) **2. ) 

eta=l . /sqrt (2 . ) *sqrt (eps+sqrt (eps**2 .+ (18 . *sig/f ) *2 . ) ) 
goto  3 
end  if 

IF  THE  POINTS  ARE  IN  DIFFERENT  MEDIA 

if (zl .le . 0 . .and. z2 . gt . 0 . )  then 

h=-zl 

z=z2 

goto  1 

else 

h=-z2 

z=zl 

end  if 


ITERATION  FOR  THE  INCIDENCE  ANGLE  BEGINS  WITH  AN  INITIAL  GUESS 
c 

1  theta=atan (d/ (h+z/eps) ) 
do  2  it=l,nit 

eta=l . /sqrt (2 . ) *sqrt (eps+sin (theta) **2 . +sqrt ( (eps- 
%  sin (theta) **2.)**2.+(18.*sig/f)**2.)  ) 
phi=asin (sin (theta) /eta) 

2  theta=atan ( (d-z*tan (phi) ) /h) 
c 

C  "nit"  th  ITERATION  VALUES  FOR  THE  INCIDENCE  AND  THE  REFRACTION  ANGLES 

C  "theta"  AND  "phi"  AND  THE  REFRACTIVE  INDEX  "eta"  ARE  TAKEN  TO  BE  THEIR 

c  SOLUTIONS. THEN  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  PARTS  OF  THE  RAY 

c  PATH  FOLLOW 

C 

dl=h/cos (theta) 
d2=z/cos (phi) 
return 
end 
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3.2  GENERAL  INSTRUCTIONS  AND  COMMENTS 


This  program  processes  the  measured  I  and  Q  data  with  the  computed  propagation 
times  to  form  the  image  of  the  specified  target  area.  It  reads  three  input  files:  First  one, 
named  image. in,  contains  parameters  and  the  names  of  the  two  other  input  files.  The 
second  one  contains  the  total  propagation  times  computed  by  the  program  proptime. f.  The 
third  one  contains  the  measured  I  and  Q  data.  It  also  makes  use  of  the  subroutine  interface 
in  calculating  the  reference  propagation  time  for  calibration  purposes. 

The  IBM  PC  version  of  this  program  is  the  PROGRAM  #4  to  be  described  next. 

3.2.1  INPUT  FILE  #1:  image.in 

This  input  file  contains  the  input  parameters  and  the  names  of  the  two  other  input  files. 


Format: 


Line# 

Input 

Parameter 

Format 

Description  of  the  Parameter 

1 

nel 

i3 

Number  of  Receiver  Array  Elements 

2 

nsamp 

ti 

Number  of  Samples  of  the  I  and  Q  waveforms 

3 

nz 

tt 

Number  of  Scan  Points  in  z  Direction 

4 

ny 

t» 

"  y  " 

5 

nx 

ti 

"  x  "  (if  no  scanning  is 

6 

iref 

i3 

desired  in  x  direction,  enter  nx=l) 

The  sample  number  corresponding  to  the  center  of  the  pulse 

7 

xtr 

fl0.4 

scattered  from  the  reference  scattering  point  and  measured  at 
the  reference  measurement  point. 
x-Coordinate  of  the  transmitter  [cm] 

8 

ytr 

it 

y 

9 

ztr 

ii 

z 

10 

xrs 

ti 

x-Coordinate  of  the  reference  scattering  point  [cm] 

11 

yrs 

it 

y 

12 

zrs 

it 

z 

13 

xrm 

ii 

x-Coordinate  of  the  reference  measurement  point  [cm] 

> 
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14 

yrm 

it 

y 

15 

zrm 

ft 

z 

16 

eps 

n 

Relative  dielectric  constant  of  the  subsurface  medium  (>1) 

17 

sig 

it 

Conductivity  of  the  subsurface  medium  [mho/meter]  (>0) 

18 

f 

n 

Frequency  of  radiation  [GHz] 

19 

tsamp 

it 

Sampling  time  interval  for  I  and  Q  data  [nsec] 

20 

iqname 

a25 

Name  of  the  input  file  containing  the  measured  I  and  Q  data 

21 

tmname 

it 

Name  of  the  input  file  containing  the  computed  propagation 

times 

22 

imname 

ti 

Name  of  the  output  file  for  storing  the  image  data 

3.2.2  INPUT  FILE  #2:  (to  be  named  by  the  user) 

This  input  file  contains  the  Inphase  (I)  and  Quadrature  (Q)  components  of  the 
measured  complex  waveforms  at  the  output  of  the  receiver  array  elements. 


Format: 

Starting  with  the  first  receiver  element,  for  each  receiver  element,  five  I&Q  pairs  per 
line  (ten  numbers  per  line),  each  number  is  a  five  digit  integer  including  its  sign.  First  the 
inphase  (I)  then  the  quadrature  (Q)  component  of  the  pair  is  stored.  There  are  500  I&Q  pair 
samples  for  each  receiver  element. 


3.2.3  INPUT  FILE  #3:  (to  be  named  by  the  user) 

This  input  file  contains  the  propagation  times  which  have  already  been  computed  by 
the  program  proptime. f  for  the  given  imaging  setup  and  the  media. 

Format: 

Format  of  the  data  can  be  easily  acquired  from  the  code  of  the  previously  discussed 
program  proptime.f.  For  each  scan  point,  as  many  propagation  times  as  the  number  of 
receiver  elements  are  stored  in  as  many  lines  of  the  input  file  as  they  can  fill.  Each  line 


contains  eight  numbers  except  for  the  last  one  which  may  contain  less.  For  the  next  scan 
point  the  storing  starts  from  the  beginning  of  the  next  blank  line.  Numbers  are  stored  in  a 
floating  number  space  with  format  "f9.5". 


3.2  DESCRIPTION  OF  IMAGE  PROCESSING 


Figure  5:  Imaging  Geometry 


Consider  Figure  5  which  shows  the  transmitter  beam  illuminating  a  strip  on  the 
ground.  Let  there  be  N  isotropic  (energy  scattered  uniformly  in  all  directions)  point 
scatterers  within  this  volume.  Each  scatterer  reradiates  a  portion  of  the  incident  wave  (i.e. 
acts  as  the  center  of  a  propagating  spherical  wave  front),  and  whose  action  is  represented 
by  an  amplitude  scattering  coefficient  and  a  time-delay  Tn  of  the  incident  pulse,  n  =  1, 2, 
....  N.  The  receiver  array  consists  of  M  elements  and  fm(t)  represents  the  received  field 
amplitude  at  the  m*  receiver  array  element.  The  transmitted  pulse  is  assumed  to  be  as  given 
in  Figure  6  where  the  envelope  is  modelled  to  be  Gaussian  (to  represent  the  band-limiting 
effect  of  the  hardware  on  the  ideal  square  pulse),  the  angular  radiation  frequency  is  to  and 
the  effective  pulse  duration  (spacing  between  3  dB  drop  points)  is  T.  Then  the  transmitted 
pulse  can  be  expressed  as 


E(t)  =  Exp 


-(In  4) 


't  ^ 


Cos  (cot) 


(7) 


Figure  6:  Transmitted  Pulse 


Further  assume  that  the  point  scatterers  cause  no  differential  envelope  delay  during 
the  process  of  scattering,  and  that  the  pulses  arrive  at  the  receiver  without  dispersion.  Then, 
using  (7)  and  the  homogeneity  of  the  medium,  we  can  express  the  received  real  signal  at  the 
output  of  the  mth  element  as 


N 

fm(0  =  ]£  Rn  ExP 
n=l 


-  (In  4) 


t-T, 


nm 


Cos[co(t  Tn-xnm)] 


,  m  =  1,  2 . M 


(8) 


where 

R,j:  Real  scattering  coefficient  of  the  n*  scatterer. 

Tnm:  Total  propagation  time  from  the  transmitter  to  the  m*  receiver  element  via  the 
n*  scatterer. 

Tn:  Amount  of  time  delay  of  the  sinusoid  of  the  pulse  upon  being  scattered  by  the 
n*  scatterer. 

co:  Carrier  angular  frequency. 


The  received  signal  fm(t)  is  a  bandpass  signal  which  can  be  represented  by 


fm(t)  =  am(t)  Cos  [cot  -  cpm(t)] 


(9) 
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4) 


which  is  a  sinusoidal  signal  with  an  angular  frequency  co,  and  time- variant  amplitude  and 
phase.  Then  let  us  define  complex  analytical  signal  as 


fm(t)  =  am(t)  exp[i<pm(t)]  (10) 

which  contains  the  same  information  as  fm(t)  does  (i.e.,  amplitude  and  phase  information.) 
Mathematically,  it  is  alright  to  write  (10)  as 


fm(0  =  am(t)  Cos  [<pm(t)]  +  i  am(t)  Sin[(pm(t)] 


(11) 


The  signal  that  we  measure  as  voltage  at  the  output  of  the  m*  receiver  antenna  is  the  real 
signal  fm(t)  represented  in  (9).  If  we  expand  the  Cosine  function  in  (9),  we  get 


fm(t)  =  |am(t)  Cos[<Pm(0]}  Cos(o)t)-{am(t)  Sin[<Pm(t)]}sin(ox) 


(12) 


Let  Im(t)  =  MO  Cos  [cpm(0]  (13a) 

and  Qm(0  =  am(t)  Sin  [<pm(0]  ( 1 3b) 

Then 

fm(t)  =  U0  Cos  (cot)  -  Q^t)  Sin  (cot)  (14) 

These  new  functions  Im(t)  and  Qm(t)  are  called  the  inphase  and  the  quadrature 
components,  respectively,  of  the  real  signal  fm(t).  From  (11)  and  (13a,b),  we  see  that  that 
Im(t)  and  Qm(t)  are  also  the  real  and  the  imaginary  parts  of  the  complex  analytical  signal 
fm(0.  That  is,  the  inphase  and  the  quadrature  components  of  a  real  bandpass  signal  are  the 
real  and  the  imaginary  parts,  respectively,  of  its  complex  signal  equivalent  which  contains 
all  the  necessary  information  about  the  real  bandpass  signal. 

Then,  in  order  to  represent  the  measured  signal  (assumed  to  be  bandpass)  in  a 
complex  signal  form,  all  we  need  to  do  is  measure  the  I  and  Q  components  of  the  signal. 
From  (14),  we  see  that  to  get  I  and  Q,  we  can  multiply  the  real  signal  with  the  reference 
signals  Cos(cot)  and  Sin(cot),  respectively,  and  then  low  pass  filter  to  get  rid  of  the  signals 
with  2co  frequency.  After  measuring  I  and  Q  components,  using  (12)  and  (13a, b),  the 
received  complex  signal  is 


fm(t)  =  Im(t)  +  i  Qm(t) ,  m=l,2 . M 


(15) 


Imaging  the  scatterers  means  estimating  the  real  scattering  coefficients  Rn>  n=  1.2....N 
by  processing  the  measured  signals  f£(t),  m=l,2,..,M.  How  the  information  about  the 
scattering  coefficients  are  contained  in  the  measurements  is  given  to  us  by  our  modeling  of 
the  received  signal  in  (8).  Since  it  is  already  in  the  form  of  (9),  its  complex  equivalent 
signal  can  be  readily  written  as 

exp[ico(Tn  +  Xjuj,)]  ,  m=l,  2,...,  M  (16) 

Define  a  complex  scattering  coefficient  Sn  for  the  n*  scatterer  as 

=  Rn  exp(icoTn) 

Then,  (16)  becomes 

exp(icoTnm)  ,  m  «  l,  2,...,  M  (17) 


fCm(t) 


>  RnExp 


(In  4) 


t-T, 


nm 


The  goal  of  the  image  reconstruction  process  is  to  obtain  "good"  estimates  of  Sn, 
n=l,2,...,N,  from  the  measurements  frmfr),  m=l,2,..,M.  A  suitable  candidate  is  the  simple 
"delay  and  sum"  beamformer  as  argued  below.  Assume  that  we  want  to  image  the  point 
(x0,y0,z0)  where  the  ngN1  poinj  scatterer  is  located,  i.e.,  we  wish  to  estimate  Sno  in  the 
presence  of  other  scatterers.  Then,  fcm(t=Tnom)  will  tend  to  have  higher  contribution  from 
Sno  comPared  to  frjrXuT^m)  as  can  be  seen  by  writing 


^m^=^n0m  )  “  ^n0exp(^^n0m)  X  Sn©xp[-(ln4)(^n0m  ^nm)  exp(itOXnm)  (18) 

n*n0 


for  m  =  1,  2,...,  M. 


A  simple  and  obvious  choice  of  an  estimate  of  Sno  from  (18),  is 


Since  we  have  access  to  the  I  and  Q  components  of  the  measurements  f£(0, 
m=i,2,..,M,  we  should  rewrite  (19)  as 


and 


sn0~ 


_1 

M 


I  [imCWn)  +  i  QmC^nom)]  exp(-  i  COtnom) 
m=l 


(20) 


There  is  one  problem  in  practice.  To  estimate  S^,  we  pick  the  best  candidates  as 

fcm(t=  Tnom),  ms  1,2 . M.  In  doing  it,  we  argued  that  these  numbers  have  higher 

contribution  from  Sno  compared  to  fcm(t  *Tnom)  since,  according  to  the  model,  the  pulse 
coming  from  the  n0th  point  scatterer  peaks  at  t=Tnom-  However,  this  is  based  on  the 
assumption  (in  the  model)  that  the  peak  of  the  pulse  leaves  the  transmitter  at  t  =  0.  One 
might  argue  that  this  is  alright  since  as  long  as  we  know  the  relative  start  time  of  the 
measurements  with  respect  to  the  firing  time  of  the  transmitter.  However,  in  practice, 
measurement  of  this  relative  time  is  not  realistic  and  practical. 


Instead,  we  have  a  more  secure  way  around  this  problem.  First  of  all,  we  keep  the 
starting  time  of  the  sampling  constant  relative  to  the  firing  of  the  transmitter  (which  is 
unknown).  We  put  a  reference  target  (a  metal  plate)  on  the  surface  and  look  at  the  amplitude 
waveform  of  the  received  signal  at  the  output  of  one  of  the  receive  array  elements.  Then, 
we  locate  in  the  time  axis  where  the  pulse  from  the  reference  target  peaks  and  call  it  tref 
("iref '  in  the  program  image. f)t  Next,  we  calculate  the  total  propagation  time  of  the  same 
pulse  from  the  transmitter  to  the  receiver  via  the  reference  target  and  call  it  Tref  ("tref'  in  the 
program).  Then  we  have  to  add  to  every  total  propagation  time  Tnm  the  amount  tref~  Tref  in 
order  to  locate  the  pulses  in  the  measured  signals  correctly.  Considering  this  correction  for 
the  calculated  propagation  times,  (20)  takes  the  final  form 


M 

I 

m=lL 


(W+'Wn 


-  Xref)  +  i  QmOref  +  ^01  ~  ^ref)]  exp  (~  i  C0Tnom) 


We  did  not  add  the  time  correction  to  the  phase  term  since  the  error  in  phase  shifting 
term  will  remain  constant  for  all  m=l,2,...,M,  and  can  be  taken  to  the  front  of  the 
summation,  and  therefore, will  not  effect  the  absolute  value  of  the  summation. 


oooooooooooooo 
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PROGRAM  #4  (IBM  PC  Version)  :  image.for 


Location:  ozdemir/research/ibm 


4.1  CODE 


c 

c 

c 

c  ********************  IBM  VERSION  **************************** 

C  ************************************************************* 

c 

c 

C  THIS  PROGRAM  FORMS  THE  IMAGE  OF  A  SPECIFIED  TARGET  REGION  USING  THE 

C  PRE-CALCULATED  PROPAGATION  TIMES  AND  THE  TIME  SAMPLES  OF  THE  INPHASE  AND 

C  THE  QUADRATURE  COMPONENTS  OF  THE  COMPLEX  WAVEFORMS  MEASURED  AT  THE 

C  OUTPUT  OF  THE  RECEIVER  ARRAY  ELEMENTS.  THE  IMAGE  FORMING  IS  PERFORMED  IN 

C  ACCORDANCE  WITH  THE  IMAGING  MODEL  PRESENTED  IN  "VFRC  PROGRESS  REPORT  NO. 

C  59"  (REFERENCE  1) . 

C 

c 

INPUT  PARAMETERS: 

nel :  NUMBER  OF  RECEIVER  ARRAY  ELEMENTS 

nsamp :  NUMBER  OF  TIME  SAMPLES  OF  THE  I  AND  Q  WAVEFORMS 
nz:  NUMBER  OF  SCAN  POINTS  IN  z-DIRECTION 
ny :  "  y  " 

nx:  "  x  " 

iref:  THE  SAMPLE  NUMBER  CORRESPONDING  TO  THE  CENTER  OF  THE  PULSE 

SCATTERED  FROM  THE  REFERENCE  SCATTERING  POINT  AND  MEASURED  AT  THE 
REFERENCE  MEASUREMENT  POINT, 
xtr:  x-COORDINATE  OF  THE  TRANSMITTER  [cm] 
ytr:  y  " 

ztr:  z  " 

xrs :  x-COORDINATE  OF  THE  REFERENCE  SCATTERING  POINT  [cm] 
c  yrs :  y  " 

c  zrs:  z  " 

c  xrm:  x-COORDINATE  OF  THE  REFERENCE  MEASUREMENT  POINT  I  cm] 

c  yrm:  y  " 

c  zrm:  z  " 

c  eps:  RELATIVE  DIELECTRIC  CONSTANT  OF  THE  SUBSURFACE  MEDIUM  (>=1) 

c  sig:  CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter]  (>=0) 

c  f:  FREQUENCY  OF  RADIATION  [GHz] 

c  tsamp:  SAMPLING  TIME  INTERVAL  FOR  I  AND  Q  DATA  [nsec] 

c  iqname :  NAME  OF  THE  INPUT  FILE  CONTAINING  THE  MEASURED  I  AND  Q  DATA 

c  tmname:  NAME  OF  THE  INPUT  FILE  CONTAINING  THE  COMPUTED  PROPAGATION  TIMES 

c  imname :  NAME  OF  THE  OUTPUT  FILE  FOR  STORING  THE  PRODUCED  IMAGE  DATA 

c 
c 

c  VARIABLES : 

c 

c  in:  INPHASE  COMPONENET  OF  THE  RECEIVED  FIELD 

c  qd:  QUADRATURE  COMPONENET  OF  THE  RECEIVED  FIELD 

c  tl:  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE  REFERENCE  SCATTERER 

c  POINT  [nsec] 

c  t2 :  PROPAGATION  TIME  FROM  THE  REFERENCE  SCATTERER  POIUT  TO  THE  REFERENCE 

c  MEASUREMENT  POINT  [nsec] 

c  trpf :  RF^FRFNOF  TOTAT.  ppopaoatton  ttmf  FROM  thp  ttfp  TO  THF, 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 
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c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


REFERENCE  MEASUREMENT  POINT  (tl  +  t2)  [nsec] 
w:  ANGULAR  FREQUENCY  OF  RADIATION  [rad/nsec] 

t:  TOTAL  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE  SPECIFIED  SCAN 

POINT  AND  THEN  TO  THE  SPECIFIED  RECEIVE  ARRAY  ELEMENT  POSITION  [nsec] 
isamp:  SAMPLE  NUMBER  CORRESPONDING  TO  THE  CENTER  OF  THE  PULSE  SCATTERED 
FROM  THE  SPECIFIED  SCAN  POINT  AND  MEASURED  AT  THE  SPECIFIED 
RECEIVE  ARRAY  ELEMENT 

s:  ESTIMATE  OF  THE  REAL  SCATTERING  COEFFICIENT  OF  THE  SCATTERER  ASSUMED 
TO  BE  AT  THE  SPECIFIED  SCAN  POINT  (VALUE  OF  THE  PRODUCED  GRAY  SCALE 
IMAGE  AT  THE  SPECIFIED  SCAN  POINT) 


dimension  t(22),s(101) 

integer  in (22, 301) ,  qd (22,  301) , in junk, qdjunk 
complex  sum 

character*25  iqname, tmname, imname 

OPEN  THE  INPUT  FILE 

open (11, file=' image . in' ) 

READ  THE  INPUT  FILE 

read(ll,'  (i3) ' )  nel, nsamp, nz , ny, nx 
read(ll, '  (i3) ')  iref 

read (11, '  (f 10 . 4 ) ' )  xtr,ytr,ztr,xrs,yrs, zrs, xrm, yrm, zrm,  eps, sig, 

%  f,tsamp 

read(ll,'  (a25)')  iqname, tmname, imname 

CALCULATE  THE  TOTAL  PROPAGATION  TIME  FROM  THE  TRANSMITTER  TO  THE 
REFERENCE  MEASUREMENT  POINT  VIA  REFERENCE  SCATTERING  POINT 

call  interface (eps,  sig,  f ,  xtr, ytr, ztr,  xrs,  yrs,  zrs,  10,  eta, dl, d2 ) 
tl= (dl+eta*d2) /30 . 

call  interface  (eps,  sig,  f,  xrs,  yrs,  zrs,  xrm,  yrm,  zrm,  10,  eta,  dl,  d2 ) 

t2=(dl+eta*d2)/30. 

tref=tl+t2 

w=8 . *atan (1 . ) *f 

OPEN  THE  INPUT  FILES  FOR  PROPAGATION  TIME  AND  IQ  DATA,  AND  THE  OUTPUT 
FILE  FOR  STORING  THE  IMAG&  DATA 

open (21, file=tmname) 
open (31, f ile=iqname) 

READ  THE  IQ  INPUT  FILE 

do  6  i-1,22 

read(31,'  (10i5)')  (.in  junk, qdjunk,  j=l,  199)  , 

%  (in(i,j),qd(i,j),j =1,301) 
open (41, f ile=imname,  status=' new' ) 

START  FORMING  THE  IMAGE 

do  3  ix=l,nx 
do  3  iz=l,nz 
do  4  iy=l,ny 

READ  THE  PROPAGATION  TIMES  THAT  TAKE  THE  WAVE  TO  TRAVZL  FROM  THE 
TRANSMITTER  TO  ALL  THE  RECEIVER  ELEMENTS  VIA  THE  SPECIFIED  SCAN  POINT 

read(21, ' (8f9.5)')  (t (iel) , iel-1, nel) 
sum=cmplx (0 . ,  0 . ) 

PHASE  SHIFTING  AND  SUMMING  OVER  RECEIVER  ELEMENTS 
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do  5  iel=l,nel 

isamp=iref +nint ( (t (iel) -tref ) /tsamp) 

^  if (isamp . gt . nsamp)  isamp=nsamp 

5  sum=sum+cmplx (float (in (iel, isamp-199) ) , 

%  float (qd (iel,  isamp-199) )) *cexp (cmplx (0 -w*t (iel) ) ) 

4  s (iy) =cabs (sum) /float (nel) 

*  c 

c  ASSIGN  TO  THE  CORNER  OF  THE  IMAGE  THE  NORMALIZATION  VALUE 
c 

if (iz .eq.nz .and.ix.eq.nx)  s (ny) =150 . 

3  write (41, ' (10li4) ' )  (nint (s (iy) ) , iy=l, ny) 
stop 
end 
c 
c 

c  THIS  SUBROUTINE  CALCULATES  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE 

C  PARTS  OF  THE  RAY  PATH  BETWEEN  TWO  SEPERATE  POINTS  (PI  AND  P2)  IN  A  MEDIA 

C  COMPOSED  OF  AIR  AND  A  SUBSURFACE  MEDIUM,  EACH  OF  WHICH  IS  A  HOMOGENOUS, 

c  ISOTROPIC,  NONMAGNETI ZABLE  AND  CONDUCTING  (EXCEPT  FOR  AIR,  WHICH  IS 

C  NONCONDUCTING)  MEDIUM,  SEPERATED  BY  A  PLANAR  INTERFACE.  THE  POINTS  CAN  BE 

c  IN  DIFFERENT  MEDIUMS  AS  WELL  AS  IN  THE  SAME  MEDIUM.  THE  INTERFACE  IS  THE 

C  Z=0  OR  (X-Y)  PLANE.  POINTS  WITH  Z<=0  ARE  IN  THE  AIR  MEDIUM  WHILE  POINTS 

c  WITH  Z>0  ARE  IN  THE  SUBSURFACE  MEDIUM, 

c 
c 

THE  PARAMETERS  AND  THE  VARIABLES  USED  IN  THE  PROGRAM: 


C 

C 

C 

c 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


eps 

(INPUT 

PARAMETER) 

sig 

( 

II 

) 

f 

( 

II 

) 

xl 

( 

II 

) 

yi 

( 

II 

) 

zl 

( 

II 

) 

x2 

( 

It 

) 

y2 

( 

II 

) 

z2 

( 

II 

) 

nit 

( 

II 

) 

eta 

(OUTPUT  PARAMETER) 

dl 

( 

II 

) 

d2 

( 

tl 

) 

RELATIVE  DIELECTRIC  CONSTANT  OF  THE  SUBSURFACE 
MEDIUM  (>=1) 

CONDUCTIVITY  OF  THE  SUBSURFACE  MEDIUM  [mho/meter] 
FREQUENCY  OF  RADIATION  [GHz] 

X-COORDINATE  OF  THE  FIRST  POINT  (PI)  [cm] 

Y  " 

Z  " 

X-COORDINATE  OF  THE  SECOND  POINT  (P2)  [cm] 

Y  •' 

Z  " 

NUMBER  OF  ITERATIONS  IN  CALCULATING  THE  INCIDENCE 
ANGLE 

REFRACTIVE  INDEX  OF  THE  SUBSURFACE  MEDIUM 
EVALUATED  AT  THE  CALCULATED  INCIDENCE  ANGLE 
THE  LENGTH  OF  THE  AIR  PART  OF  THE  RAY  PATH  [cm] 
THE  LENGTH  OF  THE  SUBSURFACE  PART  OF  THE  RAY  PATH 


[cm] 

theta  (INTERMEDIATE  VARIABLE)  :  INCIDENCE  ANGLE  [rad] 
phi  (  "  )  :  REFRACTION  ANGLE  [rad] 

d  (INTERMEDIATE  PARAMETER) :  HORIZONTAL  DISTANCE  BETWEEN  THE  POINTS  PI 

AND  P2  [cm] 

h  (  "  )  :  HIGHT  OF  THE  POINT  (PI  OR  P2)  IN  THE  AIR 

FROM  THE  INTERFACE  [cm] 

2  (  "  ) :  DEPTH  OF  THE  POINT  IN  THE  SUBSURFACE 

MEDIUM  FROM  THE  SURFACE  (INTERFACE)  [cm] 


subroutine  interface (eps, sig, f, xl,yl, zl,x2,  y2,  z2,  nit,  eta,  dl,d2) 
real  eps, sig, f , xl,yl,  zl, x2, y2, z2, eta, dl, d2 
integer  nit 

REDUCTION  OF  THE  3-D  PROBLEM  TO  A  PLANAR  GEOMETRY  PROBLEM 


d=sqrt ( (xl-x2)**2.+(yl-y2)**2.) 

IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  AIR  MEDIUM  OR  e?s=l 

if ( (zl . le . 0 . .and. z2 . le . 0 . ) .or .eps .eq. 1 . )  then 
dl=scrrt  (d**2.  +  (zl-z2)**2.) 
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d2=0  . 
eta=l . 

>  '  g°to  3 

end  if 

c 

c  IF  BOTH  POINTS  (PI  AND  P2)  ARE  IN  THE  SUBSURFACE  MEDIUM... 

*  c 

if  (zl .ge . 0 . .and. z2 .ge . 0 . )  then 
dl=0 . 

d2=sqrt (d**2 .+ (zl-z2) **2.) 

eta=l . /sqrt (2 . ) *sqrt (eps+sqrt (eps**2 .+ (18 . *sig/f ) *2 . ) ) 
goto  3 
end  if 
c 

c  IF  THE  POINTS  ARE  IN  DIFFERENT  MEDIA 

c 

if  (zl . le . 0 . .and. z2 .gt . 0 . )  then 
h=-zl 
z=z2 
goto  1 
else 
h=-z2 
z=zl 
end  if 
c 

c  ITERATION  FOR  THE  INCIDENCE  ANGLE  BEGINS  WITH  AN  INITIAL  GUESS 

C 

1  theta=atan (d/ (h+z/eps) ) 
do  2  it=l,nit 

eta=l . /sqrt (2 . ) *sqrt (eps+sin (theta) **2 .+sqrt ( (eps- 
%  sin (theta) **2.)**2.+(l8. *sig/f ) **2 . ) ) 
phi^asin (sin (theta) /eta) 

2  theta*=atan  ( (d-z*tan  (phi) ) /h) 
c 

C  "nit"  th  ITERATION  VALUES  FOR  THE  INCIDENCE  AND  THE  REFRACTION  ANGLES 

C  "theta"  AND  "phi"  AND  THE  REFRACTIVE  INDEX  "eta"  ARE  TAKEN  TO  BE  THEIR 

C  SOLUTIONS. THEN  THE  LENGTHS  OF  THE  AIR  AND  THE  SUBSURFACE  PARTS  OF  THE  RAY 

C  PATH  FOLLOW 

c 

dl=h/cos (theta) 
d2=z/cos  (phi)  «•' 

return 
end 
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4.2  GENERAL  INSTRUCTIONS  AND  COMMENTS 


This  program  is  the  modified  version  of  PROGRAM  #3  to  meet  the  memory 
allocation  requirements  of  IBM  PCs.  Every  thing  else  is  exactly  the  same  as  for 
PROGRAM  #3  (input,  output  files,  formats,  etc.) 

The  runing  time  of  the  program  is  roughly  proportional  to  the  number  of  scan 
points  times  the  number  of  receiver  elements.  Let  the  number  of  scan  points  in  longitudinal, 
vertical  and  horizontal  directions  be  Nx,  Nz,  Ny,  and  the  number  of  receiver  elements  be 
Nr,  respectively.  When  Nx=l,  Nz=51,  Ny=101  and  Nr=22,  the  program  took  14.5  minutes 
to  run.  Then  the  formula  for  runing  time  (on  IBM  PC)  is 

Runing  Time  *  14.5  x  (  Nx  x  Nz  x  Ny  x  Nr )  /  ( 1  x  51  x  101  x  22  )  Minutes. 
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(Orientation  of  Metal  Mine  Casing) 


Figure  2  :  Isometric  Image  of  Plastic  Mine  3”  (7.5  cm)  Below  the  Surface 
Er=  3.2,  o  =  0.005  Date:  12/06/89 
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Figure  4  :  Isometric  Image  of  Plastic  Mine  3"  (7.5  cm)  Below  the  Surface 
Ej  =  5.5,  a  =  0.145  Dale:  08/01/90 


